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20.  ABSTRACT 

mined.  The  simulation  of  the  two-stage  fan  component  was  carried  out  using 
a quasi  one-dimensional,  pitch  line,  blade-row  by  blade-row  model  which 
used  blade-row  loss  coefficients  and  deviation  angles  determined  from 
analysis  of  clean-inlet-flow  test  data.  This  model  was  ideally  suited  for 
studying  the  fan  response  to  the  planar  waves  produced  by  the  Planar  Pressure 
Pulse  Generator,  a device  which  produced  highly  planar  waves  with  peak-to- 
peak  fluctuations  of  10-30  percent  of  the  mean  total  pressure  over  a frequency 
(range  of  4?  to  800  Hz. 

The  results  of  the  simulations  showed  that  the  model  accurately  reproduced 
(the  waveshape,  amplitude,  and  phase  of  the  experimentally  determined  fan 
'Inlet  and  exit  pressures  and  the  fan  inlet  flow  up  to  75  Hz  where  the  need 
for  including  unsteady  blade  circulation  effects  on  blade  characteristics 
began  to  become  apparent.  Further,  the  blade  row  model  showed  accurate 
simulation  of  the  pressures  in  the  interconnecting  ducts  by  reproducing  the 
wave  and  its  harmonics  up  to  approximately  118  Hz  where  the  length  of  the 
volumes  becomes  too  large  to  provide  proper  spatial  definition  of  the 
harmonics  of  the  standing  waves  although  the  fundamental  frequency  is 
properly  reproduced.  In  addition,  the  model  accurately  predicted  the 
trends  and  approximately  predicts  the  level  of  fan  dynamic  distortion 
sensitivity  and  transfer. 

A method  for  including  the  unsteady  blade  circulation  effects  in  the 
fan  simulation  has  been  investigated  and  appears  to  improve  the  phase 
relationships  between  predicted  and  test  values  at  higher  frequencies. 

Further  refinement  of  this  method  is  required. 
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1 . 0 SUMMARY 


The  objectives  of  this  program  were  to  demonstrate  that  a one-dimensional 
pitch-line  Dynamic  Digital  Blade  Row  Compression  Component  Stability  Model 
could  provide  an  accurate  simulation  of  a fan  component  when  the  fan  was 
subjected  to  planar  wave  oscillations  and  to  show  the  manner  in  which  the 
fan  aerodynamics  responded  to  the  unsteady  flows. 

The  capability  of  the  model  to  provide  accurate  simulations  was  demon- 
strated by  conducting  simulations  of  test  points  obtained  during  testing  of 
a two-stage  fan  component  and  the  Planar  Pressure  Pulse  Generator  used  to 
generate  an  unsteady-planar-wave  fan  inflow.  The  capability  of  the  model 
to  handle  unsteady  flows  was  improved  by  introducing  an  analytical  expres- 
sion for  the  mean  pressure  acting  on  the  lateral  surfaces  of  a blade  free 
volume,  by  including  the  tangential  force  resulting  from  the  storage  of 
angular  momentum  within  a blade  row  volume,  and  by  developing  a method  for 
handling  the  unsteady  effects  upon  blade  circulation. 


The  two-stage  fan  map  was  accurately  simulated  at  80%  and  100%  corrected 
speeds.  Accurate  prediction  of  the  surge  line  was  made  through  use  of  a 
previously  established  technique  which  is  based  upon  the  ratio  of  the  flow 
leaving  a blade  row  volume  to  the  flow  entering  the  volume. 

Simulations  of  the  Planar  Pressure  Pulse  Generator  test  data  were  made 
over  the  42  Hz  to  350  Hz  frequency  range.  These  simulations  gave  accurate 
replication  of  the  measured  waveforms  throughout  the  model  over  the  42  Hz 
to  80  Hz  frequency  range  at  both  80-  and  100-percent  corrected  speeds.  At 
higher  frequencies  the  model  had  increasing  difficulty  in  duplicating  the 
harmonics  in  the  waveforms  such  that  at  220  Hz,  the  waveform  harmonics  are 
completely  washed  out,  although  the  fundamental  frequency  appears  to  be 
adequately  reproduced.  The  lack  of  reproduction  of  the  harmonics  is  attri- 
buted to  the  presence  of  standing  waves  and  the  lack  of  spatial  resolution; 
that  is,  the  selected  volume  lengths  are  not  sufficiently  small  compared  to 
the  wavelengths  associated  with  the  harmonic  frequencies  that  are  present. 
Planar  wave  distortion  transfer  was  calculated  from  the  model  results  and 
its  value  was  found  to  fall  between  the  hub  and  tip  test  values  as  should 
be  the  case  for  a pitch  line  model.  In  fact,  the  model  predicted  the 
anomolous  amplification  that  was  present  in  the  test  data  at  118  Hz,  80% 
corrected  speed. 

Analysis  of  the  stationary  data  prior  to  throttling  to  surge  showed 
that  at  frequencies  greater  than  80  Hz  (w*  s 0.25),  the  model-calculated- 
parameters  downstream  of  the  fan  were  out  of  phase  with  the  measured  vari- 
ables with  the  former  leading.  The  phase  angle  also  appeared  to  be  a func- 
tion of  frequency.  This  behavor  was  anticipated  since  the  basic  model  for- 
mulation assumed  that  the  blade  characteristics  adjusted  instantaneously  to 
the  new  incidence  angle.  Inclusion  of  an  effective  incidence  angle  based 
upon  a blade-row  time  constant  led  to  results  which  produced  a lag  in  the 
calculated  downstream  parameters.  However,  further  refinement  of  the  time 
constant  is  needed. 
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Prior  to  discussing  the  results  obtained  by  throttling  the  model  to 
surge,  it  should  be  noted  that  the  surge  line  results  are  based  upon  blade 
characteristics  which  have  been  extrapolated  beyond  values  associated  with 
the  steady-state  surge  line.  Such  extrapolations  are  necessary  if  the 
model  is  to  function  in  the  presence  of  unsteady  flows.  The  model  was 
throttled  to  surge  with  the  planar-wave  inlet  conditions  and  several  interest- 
ing f indings  emerged.  It  was  found  that  the  planar-wave  distortion-sensitivity 
trends  and  magnitudes  were  essentially  reproduced  out  to  75  Hz,  but  at  118 
Hz  and  above,  especially,  for  100%  corrected  speed,  the  magnitude  of  the 
sensitivity  deviated  considerably  from  test  values.  Further,  it  was  interest- 
ing to  note  that  at  the  experimentally  determined  surge  points,  the  model 
exhibited  behavior  that  appeared  to  duplicate  the  actual  fan  response;  that 
is,  at  100%  corrected  speed,  the  model  instability  developed  extremely 
rapidly  and  computation  ceased  abruptly  while  at  the  80%  corrected  speed, 
the  model  exhibited  a much  slower  developing  instability  marked  by  erratic 
behavior  of  the  fan  internal  aerodynamics.  One  parameter  of  particular 
interest  was  the  tangential  blade  force  associated  with  the  storage  of 
angular  momentum  within  a blade-row  volume.  This  parameter  was  very  sensi- 
tive to  developing  instabilities  and  its  increasingly  time  dependent  and 
erratic  behavior  was  the  first  evidence  that  an  aerodynamic  instability  was 
developing,  although  all  basic  flow  parameters  were  well  behaved.  It  is 
possible  that  this  force  could  provide  an  important  analytic  link  between 
the  aerodynamic  and  aeroelastic  responses  of  a compression  component. 

The  results  of  this  study  are  most  encouraging  with  respect  to  demon- 
strating the  validity  of  the  Dynamic  Digital  Blade  Row  Compression  Component 
Stability  Model.  Improvements  to  this  model  can  be  made  in  three  areas: 

1)  Volume-length  selection  with  respect  to  increasing  the  frequency  response 
of  the  model,  2)  Improvement  of  the  post-surge  extrapolated  blade-row 
characteristics  and,  3)  Refinement  of  the  technique  for  including  the 
finite  time  needed  for  the  blade  circulations  to  readjust. 


2.0  INTRODUCTION 


The  time-varying  characteristics  of  total-pressure  distortion  lead  to 
documented  losses  in  surge  pressure  ratio  (References  1 and  2).  In  these 
studies,  the  loss  in  surge  pressure  ratio  can  be  viewed  as  resulting  from 
an  imposed  steady-state  distortion  and  a random  time-varying  distortion 
fluctuating  about  the  steady-state  level.  While  these  studies  provide 
gross  information  concerning  the  effect  of  the  flow  unsteadiness  charac- 
teristic of  that  produced  by  aircraft  inlets,  its  random  nature  precludes 
deterministic  analysis  of  the  unsteady  aerodynamic  processes  which  take 
place  within  a compression  component  and  the  manner  in  which  they  induce 
surge.  Compression  component  response  to  total-pressure  fluctuations  as  a 
function  of  frequency  was  determined  by  subjecting  a two-stage  fan  to  inlet 
planar-wave  total-pressure  fluctuations,  that  is,  deterministic  time- 
varying  oscillations.  This  effort  has  been  reported  upon  by  Reynolds, 

Vier,  and  Collins  (Reference  3). 

Experimental  aerodynamic  stability  investigations  such  as  those  reported 
in  References  1-3  involve  considerable  expense  considering  the  costs  asso- 
ciated with  fabricating  an  unsteady-flow  generator,  installing  a test 
vehicle  in  a cell,  fully  instrumenting  a test  vehicle,  and  running  the  test 
vehicle  over  a wide  range  of  test  conditions.  Even  if  there  are  no  facility 
limitations,  the  investigator  is  faced  with  instrumentation  accuracy  limita- 
tions which  often  make  proper  data  interpretation  a difficult  or  impossible 
task.  Clearly,  it  would  be  desirable  if  an  approach  could  be  found  which 
would  supplant  or  at  least  complement  such  testing,  and  remain  an  economi- 
cally viable  alternative. 

One  alternative  that  is  being  developed  by  the  Aircraft  Engine  Group 
of  the  General  Electric  Company  is  called  the  Dynamic  Digital  Blade  Row 
Compression  Component  Stability  Model.  This  procedure  models  the  compres- 
sion component  on  a quasi  one-dimensional,  pitch  line,  blade-row  by  blade- 
row,  time  dependent  basis.  The  capabilities  of  this  model  offer  an  unique 
approach  for  studying  the  effect  of  planar  wave  flows  upon  the  aerodynamic 
stability  of  compression  system  components. 

In  quoting  from  another  work  (Reference  4)  "....  if  such  models  are  to 
achieve  their  ultimate  capability  as  a design  and  evaluation  tool,  it  is 
necessary  to  develop  and  validate  the  capability  of  the  models  against 
existing  bodies  of  test  data.  This  step  establishes  confidence  in  the 
model,  helps  to  define  its  range  of  validity,  and  insures  that  it  is  under- 
stood." The  impact  of  this  statement  cannot  be  underestimated.  Although  a 
number  of  time-dependent  computer  models  have  been  developed  (References  5, 

6,  and  7),  the  investigators  either  did  not  possess  or  have  access  to 
large-amplitude,  experimental  data  over  a wide  frequency  range.  Such  data 
include  the  nonlinear  aerodynamic  interactions  which  occur  within  a compres- 
sion component  and  provide  the  proper  benchmark  against  which  a dynamic 
model  should  be  measured.  General  Electric  acquired  such  data  for  a two- 
stage  fan  (Reference  3)  under  United  States  Air  Force  Aero  Propulsion 
Laboratory  Contract  F33615-73-C-2001 . 
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With  these  data  in  hand,  Ceneral  Electric  was  able  to  embark  upon  the 
effort  documented  in  this  report.  The  two  main  objectives  to  this  program 
were:  l)  To  demonstrate  that  a one-dimensional,  dynamic  digital  model, 

consisting  of  a Planar  Pressure  Pulse  Generator  (P^G)  model  coupled  to  a 
two-stage  fan  model,  would  provide  an  accurate  simulation  of  the  previously 
acquired  test  data,  thereby  validating  the  dynamic  modeling  technique 
currently  being  utilized  by  the  Aircraft  Engine  Group  of  General  Electric 
and  to  determine  the  upper  frequency  limit  of  the  model;  and  2)  To  show  the 
manner  in  which  the  aerodynamics  of  the  fan  responded  to  the  unsteady  flow 
created  by  the  P^G  through  detailed  study  of  surge  pressure  ratio  loss 
sensitivity,  blade  row  and  overall  amplitude  transmission  coefficients,  and 
other  calculated  parameters. 

The  objectives  were  accomplished  by  simulating  the  test  conditions, 
analyzing  the  results,  and  by  doing  additional  detailed  analysis  of  some  of 
the  test  data.  Specifically,  the  one-dimensional  digital  P^G  and  two-stage 
fan  models  were  coupled  so  that  the  critical  P^G  exit/fan-inlet-flow  boundary 
conditions  were  self-established.  The  geometry  (lengths  and  areas)  of  the 
entrance  to,  of  the  connecting  and  exit  ducting  of  the  P^G,  and  of  the  two- 
stage  fan  were  included  in  the  model;  that  is,  the  test  setup  from  the 
bellmouth  entrance  to  the  bypass  discharge  valve  was  simulated.  The  boundary 
conditions  employed  during  this  study  assumed  a constant  total  pressure  and 
constant  total  temperature  at  the  bellmouth  entrance  and  a constant  static 
pressure  at  the  plane  of  the  bypass  discharge  valve. 

During  the  first  phase  of  this  program,  the  80 % and  100%  corrected 
speed  lines  for  the  two-stage  fan  were  simulated  using  steady,  clean  inlet 
flow  conditions.  This  step  was  accomplished  using  the  pitchline  relative 
total-pressure  loss  coefficients  and  deviation  angles  for  the  rotors  that  were 
previously  developed  from  clean-inlet-flow  fan  component  data  and  a General 
Electric  stage-stacking  procedure.  The  clean  surge  line  was  established 
based  upon  the  development  of  internally  generated  unsteady  flows  within 
the  fan  which  lead  to  aerodynamic  instability  and  was  shown  to  accurately 
reproduce  the  experimentally  determined  surge  line.  Upon  completion  of 
this  phase,  the  P^G  amplitudes  and  frequencies  generated  for  surge  points 
obtained  during  the  February  1973  test  of  the  P^G  and  the  two-stage  fan 
were  simulated  (i.e.,  a baseline  point  which  was  a near-surge  point  and  the 
ensuing  throttling  to  surge).  The  validity  of  the  model  was  demonstrated 
by  comparing  the  amplitude,  waveform  shape,  and  phase  of  the  calculated 
parameters  at  the  fan  entrance  and  exit.  Additional  checks  upon  the  digital 
model  are  provided  by  pressures  measured  upstream  of  the  P^G,  between  the 
P^G  and  fan  entrance,  and  in  the  discharge  system.  Simulation  of  the  pres- 
sures at  these  locations  permitted  determining  whether  a parameter  was 
acting  predominately  as  a standing  wave  or  a traveling  wave  in  the  duct 
connecting  the  P^G  to  the  two-stage  fan. 


It  was  anticipated  that  comparison  of  some  model  results  with  test  data 
would  prove  difficult  if  carried  out  directly.  This  was  due  to  the  presence 
in  the  test  data  of  low-level  electrical  noise,  aerodynamic  noise  (turbulence), 
and  higher  order  harmonics  of  the  impressed  fundamental.  The  basic  character 
of  the  waveforms  was  enhanced  over  the  originally  reported  results  (Reference 
3)  by  filtering  out  high  frequency  content  and  employing  a digital  analysis 
technique  - ensemble  averaging  - which  essentially  removed  any  random  con- 
tributions to  the  waveform.  This  technique  assumed  that  the  data  was 
stat istically  stationary.  Favorable  amplitude  and  phase  angle  comparisons 
of  the  model-produced  signatures  with  the  experimentally  determined  test 
signatures  at  each  frequency  and  corrected  speed  combination  provided  a 
strong  basis  for  considering  the  model  validated.  However,  it  was  antici- 
pated that  the  degree  of  simulation  would  necessarily  begin  to  deteriorate  for 
frequencies  in  the  range  of  100-200  Hz  based  upon  analytical  cascade  results 
(Reference  8).  This  would  be  an  indication  that  the  loss  coefficients  and 
deviation  angles  drived  from  steady-state  data  could  not  be  used  to  predict 
the  "high-frequency"  dynamic  performance  of  the  fan.  This  knowledge,  for 
the  first  time,  would  give  a realistic  value  on  the  upper  limit  of  the 
frequency  range  for  which  these  types  of  dynamic  models  are  valid. 

With  the  model  essentially  validated,  each  of  the  simulation  cases  was 
then  dissected  to  provide  detailed  information  as  a function  of  frequency 


• Overall  fan  dynamic  transfer  characteristics 

• Stage  transmission  characteristics 

• Overall  fan  surge  pressure  ratio  loss  per  unit  amplitude 

• Range  of  incidence  angles  that  each  rotor  encountered 

• Estimation  of  the  stage  where  instability  initiated 

Additionally,  effort  was  expended  towards  developing  a method  for 
including  the  effect  of  the  frequency  dependency  of  the  lift  force,  loss 
coefficients,  deviation  angles,  etc. 


This  report  provides  a detailed  account  of  the  manner  in  which  these 
objectives  were  approached  and  the  results  that  were  obtained. 


3.0  PLANAR  PRESSURE  PULSE  GENERATOR  AND 
TWO- STAGE  FAN  TEST  AND  RESULTS 


Because  the  information  derived  from  USAF  Aero  Propulsion  Laboratory 
Contract  F33615-73-C-2001  (Reference  3)  is  such  an  integral  part  of  the 
present  effort,  the  test  layout,  instrumentation,  and  results  are  reviewed  in 
this  section  to  help  provide  perspective  to  the  approach  followed  in  this 
study  and  to  obviate  the  need  to  refer  to  Reference  3 unless  more  specific 
details  are  required. 

3. L TEST  DESCRIPTION 

The  fan  component  tested  during  February  of  1973  in  the  Aircraft  Engine 
Croup's  Large  Fan  Test  Facility  located  in  Lynn,  Massachusetts  was  a two- 
stage  fan  of  flight  type  design.  This  fan  has  a hub-to-tip  radius  ratio  of 
0.486  and  a flapped  1GV  (Inlet  Guide  Vane)  trailing  edge  which  is  positioned 
as  a function  of  corrected  speed.  The  fan  flowpath  is  defined  in  Section 
4.4.1.  A clean  inlet  flow  map  for  this  fan  is  given  in  Figure  1. 

The  General  Electric-developed  Planar  Pressure  Pulse  Generator  (P^G)  was 
installed  approximately  one  fan-face  diameter  upstream  of  the  fan  face  and 
was  used  to  produce  the  large  amplitude  planar  waves  used  for  evaluating  the 
response  of  the  fan  to  deterministic  unsteady  flows.  The  p3g  is  a choked-  , 
flow  device  which  uses  a single-stage  rotor  and  stator  combination  to  sinu- 
soidally modulate  the  minimum  area.  The  frequency  of  the  planar  waves  is 
governed  by  the  rotor-to-stator  spacing.  The  advantage  of  this  device  lies 
in  its  ability  to  produce  highly  planar  waves  with  peak-to-peak  amplitudes  of 
10-30%  of  the  mean  total-pressure  over  the  420-800  Hz  frequency  range. 

A scaled  installation  layout  is  given  in  Figure  2 and  shows  the  relation- 
ship of  the  bellmouth,  the  p3g,  the  two-stage  fan,  the  exit  ducting,  and  the 
discnarge  valves  to  each  other. 

This  fan  was  highly  instrumented  since  the  planar  wave  testing  followed 
a performance  test  of  the  fan.  In  addition  to  the  standard  extensive  array 
of  steady-state  instrumentation,  the  installation  was  instrumented  with  high- 
response  wall-static  pressure  transducers,  total-pressure  probes,  and  a flow 
(pV)  probe.  The  locations  of  the  instrumentation  which  provided  useful  data 
of  interest  to  this  study  are  also  shown  on  the  installation  layout  of  Figure 
2.  The  letters  PS  contained  within  a probe  identifier  designated  wall- 
static-pressure  transducers  located  upstream  of  the  p3g,  between  the  P^G  and 
the  fan  IGV's,  and  in  the  fan  discharge  duct.  Total-pressure  transducers 
were  designated  by  PT  and  were  located  at  the  IGV's  as  was  the  flow  probe 
designated  by  WCA.  The  outputs  of  all  high-r«>cponse  instrumentation  were 
recorded  on  analog  magnetic  tape  for  permanency  and  subsequent  analyses. 


Figure 


Scaled  Installation  Layout  and  Instrumentation  Locat 


The  procedure  that  was  followed  In  operating  the  Planar  Pressure  Pulse 
Generator  and  fan  system  was  to  set  the  P3c  rotor-to-stator  spacing  at  its 
maximum  value.  The  fan  was  brought  to  its  proper  speed  and  the  pressure 
ratio  for  the  baseline  point  was  established  by  setting  the  bypass  discharge 
valves.  The  P^G  was  then  brought  to  the  desired  frequency  and  the  desired 
planar  wave  amplitude  was  set  by  closing  down  the  rotor-to-stator  spacing. 
Baseline  data  were  taken.  The  fan  was  then  throttled  to  surge  by  slowly 
closing  the  bypass  discharge  valves;  stationary  data  points  were  taken  along 
the  speed  line  in  order  to  establish  the  influence  of  planar  waves  upon  the 
pressure  ratio  flow  relationship.  Once  the  fan  surged,  the  bypass  discharge 
valves  were  opened  to  relieve  the  fan  and  to  permit  it  to  recover. 

During  this  test,  the  dynamic  characteristics  of  the  fan  were  investi- 
gated at  80-  and  100-percent  corrected  speeds.  Twenty-four  surges  were 
induced  as  follows: 

10  unsteady  flow  surges  at  80%  N//0 

13  unsteady  flow  surges  at  100%  N//e 

1 steady  clean  inlet  flow  posttest  surge  at  100%  N //9 

The  posttest  surge  was  used  as  a check  on  fan  deterioration.  Of  particular 
interest  to  the  current  study  effort,  are  the  sixteen  surges  given  in  Table  1 
along  with  the  baseline  points  that  were  to  be  used  for  waveform  fidelity 
comparisons.  The  difference  between  the  16  surges  with  unsteady  inlet  flow 
conditions  of  Table  1 and  the  23  surges  obtained  during  the  test  represents 
repeated  points.  However,  as  will  be  noted  by  examining  the  100%  N//0 
points,  there  are  wo  repeated  points,  one  at  118  Hz  and,  the  other  500  Hz. 
The  purpose  of  this  ostensible  repetition  is  to  determine  if  there  is  an 
effect  of  wave  amplitude  upon  dynamic  distortion  sensitivity. 

3.2  RESULTS  OK  TEST  DATA  ANALYSIS 

The  results  which  were  obtained  during  the  Planar  Pressure  Pulse 
Gennerator /two-stage  fan  test  and  which  are  particularly  pertinent  to  the 
current  study  are  discussed  and  reproduced  in  this  section.  The  sensitivity 
of  the  fan  to  planar  wave  fluctuations  is  given  in  Figures  3 and  4 for  80% 
and  100%  corrected  speeds,  respectively.  The  dynamic  distortion  sensitivity 
is  defined  as  the  loss  in  surge  pressure  ratio  at  constant  flow  measured 
between  the  clean  surge  line  and  the  mean  operating  point  with  a planar  wave 
amplitude  of  [ (P’l’MAX  “ pTMIn) /2Ptavg  1 • Hence,  the  dynamic  distortion  sensi- 
tivity is  the  loss  of  surge  pressure  ratio  per  unit  amplitude  of  the  planar 
wave.  This  definition  is  in  keeping  with  steady-state  distortion  sensitivity 
coefficients  which  are  defined  as  the  loss  in  surge  pressure  ratio  per  unit 
distortion  (circumferential  or  radial).  The  data  show  that  for  frequencies 
up  to  approximately  that  equal  to  fan  rotational  speed,  a loss  in  surge  pres- 
sure ratio  is  produced.  Above  that  frequency  essentially  no  loss  in  surge 
pressure  ratio  is  noted. 
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Table  1.  Summary  of  Surge  Points. 


f 

Rotor  to 
Stator  Spac ing 

Surge 

Number 

Baseline 
Reading  Number 

80 

42 

0. 100 

19 

82 

80 

75 

0.150 

18 

77 

80 

118 

0.090 

11 

61 

80 

240 

0.090 

9 

57 

80 

350 

0.090 

7 

53 

80 

641 

0.090 

22 

91 

80 

800 

0.060 

23 

92 

100 

42 

0.150 

16 

71 

100 

80 

0.150 

14 

69 

100 

118 

0.150 

12 

64 

100 

118 

0.090 

13 

64 

100 

220 

0.150 

2 

41 

100 

350 

0.150 

3 

43 

100 

500 

0.150 

4 

46 

100 

500 

0.091 

5 

47 

100 

613 

0.090 

20 

86 

I<> 


The  manner  in  which  the  fan  amplifies  or  attenuates  the  planar  waves  is 
described  by  the  planar  wave  distortion  transfer  coefficient.  This  coef- 
ficient is  the  ratio  of  the  nondimens ionalized  peak-to-peak  total-pressure 
fluctuations  at  the  fan  exit  divided  by  the  nond imensional ized  peak-to-peak 
total-pressure  fluctuations  at  the  fan  inlet.  Because  the  blade  loadings  are 
different  at  the  hub  and  the  tip,  it  is  expected  that  the  planar  wave  distortion 
transfer  coefficient  should  show  differences  between  values  at  the  hub  and  tip. 
The  hub  transfer  coefficient  is  shown  in  Figure  5 and  the  tip  transfer  coef- 
ficient is  shown  in  Figure  6.  It  is  interesting  to  note  that,  in  general, 
both  the  hub  and  the  tip  of  the  fan  attenuate  planar  waves  with  the  hub 
showing  pronounced  attenuation.  This  is  a beneficial  effect  for  any  com- 
pressor located  downstream  of  the  fan.  However,  it  is  worth  noting  that  both 
the  hub  and  the  tip  of  the  fan  show  an  apparent  amplification  at  118  Hz. 

Although  these  points  appeared  to  be  out  of  place  at  the  time  the  data 
analyses  were  performed,  no  justification  could  be  found  for  deleting  these 
points.  Discussion  later  in  this  report  will  return  to  these  data. 

As  a check  on  fan  "health,"  a posttest  clean-inlet-flow  surge  line 
determination  check  was  performed  at  100%  N//g".  This  check  showed  that  due 
to  deterioration,  the  fan  lost  approximately  1.2  percent  in  flow  at  surge  and 
0.017  in  surge  pressure  ratio  as  shown  in  Figure  7 during  the  planar  wave 
testing  of  the  fan. 


3 . 3 DATA  ENHANCEMENT^  ANALYSIS 

The  results  of  the  previous  section  were  in  part  obtained  by  analysis  of 
the  analog  data  recorded  on  magnetic  tape.  The  presence  of  electrical  noise, 
turbulence,  and  other  spurious  disturbances  which  are  random  tend  to  distort 
the  basic  character  of  the  waveform.  Further,  the  higher  harmonics  associated 
with  wave  reflections  within  the  test  facility  also  distort  the  data.  These 
disturbances  do  not  compromise  the  validity  of  the  results  previously  reported 
(Reference  3)  or  those  reproduced  in  Section  3.2.  However,  since  one  of  the 
primary  objectives  of  this  study  is  to  validate  the  modeling  technique  by 
comparing  computed  waveforms  with  the  test  waveforms,  it  is  necessary  that  the 
comparisons  be  carried  out  using  comparahle  bases.  The  values  calculated 
by  the  model  are  necessarily  free  from  any  random  disturbances.  Therefore, 
comparisons  using  these  results  should  be  made  with  test  data  which  also  are 
free  of  such  disturbances.  The  operation  upon  test  data  to  help  reveal  its 
basic  character  free  from  random  or  spurious  disturbances  is  termed  "Data 
Enhancement . " 

Detiiled  analysis  of  the  instantaneous  calculated  airflow  was  undertaken 
since  its  behavior  appeared  to  be  most  prejudiced  by  disturbances.  A hard 
wired  ana  log-to-d ig i t a 1 computer  (Time  Data  TDA-53P)  was  employed  for  digi- 
tizing and  enhancing  the  airflow  data  by  "ensemble  averaging."  The  instan- 
taneous airflow  waveform  was  digitized  over  several  cycles  and  each  cycle 
overlaid  on  top  of  the  next  in  such  a manner  that  an  average  waveform  is 
computed  by  the  stmimat ion  of  all  samples.  Proper  scaling  was  maintained  by 
dividing  by  the  nuriber  of  samples.  Thus,  any  disturbances  which  are  random 
become  smoothed  oi:l  essentially  eliminated  from  the  resultant  waveform  when  a 
sufficiently  large  number  of  data  samples  are  taken. 
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Figure  5*  Fan  Hub  Planar  Wave  Distortion  Transfer 


TIP  PLANAR  WAVE  DISTORTION  TRANSFER  COEFFICIENT 


Figure  8 shows  a comparison  of  the  raw  data  waveform  for  the  calculated 
airflow  in  comparison  to  the  enhanced  data  waveform.  The  ensemble  averaged 
data  indicate  that  the  fundamental  frequency  and  its  first  harmonic  are  the 
only  significant  contributors  to  the  overall  waveform  and  Lhat  all  other 
contributors  are  either  electronic  noise  or  random  aerodynamic  events.  This 
procedure  verified  that  the  primary  events  of  interest  were  confined  to  the 
first  two  or  three  harmonics  of  the  recorded  data.  Fluctuations  occurring  at 
frequencies  above  this  limit  did  not  contribute  to  the  repetitive  character 
of  the  waveform  and  could  therefore  be  eliminated  from  the  data  analysis  and 
comparisons.  As  a result  of  this  finding  and  in  view  of  the  expense  of  this 
type  of  enhancement  analysis,  an  alternate  procedure  was  sought  to  permit 
enhancement  of  all  variables  for  which  comparisons  between  test  and  model 
results  were  to  be  made.  It  was  determined  that  all  data  would  be  filtered 
at  approximately  3 times  its  fundamental  frequency  thus  eliminating  the 
spurious  noise  contributions  at  the  higher  frequencies.  Figure  8 also  shows 
the  result  of  the  filtering  procedure  in  comparison  to  the  test  data  and 
enhancement  results.  Small-amplitude,  random  disturbances  distort  the  wave 
shape  locally,  however,  the  basic  nature  of  the  waveform  is  clearly  demon- 
strated. Therefore,  the  procedure  adopted  for  data  enhancement  of  all  param- 
eters used  in  the  comparisons  was  to  filter  the  raw  data  at  a frequency  commen- 
surate with  the  third  harmonic  of  the  imposed  fundamental  P^G  frequency. 

As  a point  of  information  it  should  be  noted  that  at  the  higher  fre- 
quency test  points  (f>500  Hz),  this  procedure  did  not  always  produce  mean- 
ingful! results  with  respect  to  the  airflow  measurement.  Spurious  signals 
from  the  experimental  airflow  probe  were  present  and  were  not  eliminated  by 
the  filtering  process.  Ensemble  averaging  of  these  data  indicated  that  there 
was  no  periodic  components  for  some  readings.  For  these  readings  it  was 
concluded  that  the  probe  was  improperly  functioning.  The  ensemble  averaged 
airflow  waveform  is  shown  in  Figure  9 relative  to  the  instantaneous  pressure 
ratio  for  42  Hz  at  100%  corrected  speed.  Similarly  enhanced  waveforms  for 
other  test  points  are  shown  in  Appendix  A. 
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Figure  9-  Ensemble  Averaged  Waveforms,  Kdg.  71.  lOO%  N/\®,  *»2  Hz. 
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4.0  DESCRIPTION  OF  DYNAMIC  DIGITAL  MODEL 


The  equations  basic  to  the  Dynamic  Digital  Blade  Row  Compression  Compo- 
nent Stability  Model  are  reproduced  in  this  section  to  give  perspective  to 
the  method  employed  in  this  investigation.  To  a large  degree,  this  review 
is  taken  directly  from  References  4 and  9.  In  addition,  refinements  to  the 
model  that  have  been  developed  in  the  course  of  this  investigation  are 
detailed  as  well  as  a description  of  the  digital  model  of  the  Planar  Pressure 
Pulse  Generator.  Further,  the  geometrical  data  and  fan  characteristics  data 
are  presented,  the  externally  imposed  boundary  conditions  are  given,  and  the 
method  of  operating  the  model  is  discussed. 


OVERVIEW  OF  MODEL  TECHNIQUES 


The  complete  set  of  non-linear  partial  differential  equations  which 
describe  the  transfer  and  storage  of  mass,  momentum,  and  energy  within  a 
fluid  are  called  the  equations  of  change  (Reference  10).  These  equations 
have  been  integrated  once  over  an  arbitrary  volume  of  the  flow  system  to 
obtain  the  macroscopic  balances  for  quasi  one-dimensional  flow  without  heat 
transfer  and  are  reproduced  below  in  the  form  in  which  they  are  used  in  the 
dynamic  compression  component  model. 
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The  subscripted  variables  on  the  right-hand  side  of  the  equations  refer  to 
quantities  at  the  inlet  and  exit  of  the  control  volume.  Variables  on  the 
left-hand  side  refer  to  volume  averaged  quantities,  i.e.,  in  generalized 
form 


(4) 


The  energy  equation  (Equation  3)  was  derived  by  combining  the  equation  of 
change  for  energy  and  one  of  the  thermodynamic  T ds  relationships. 

This  set  of  equations  (other  than  being  applicable  to  quasi  one-dimen- 
sional flows  without  heat  transfer  and  to  a finite,  but  small  volume)  properly 
and  exactly  describes  the  state  of  a fluid  in  motion.  In  order  for  a solution 
to  Equations  1 to  3 be  obtained,  it  is  necessary  to  supply  the  caloric  and 
thermal  equations  of  state  and  expressions  for  Fg  (Blade  Force),  Pm  (Mean 
Pressure),  and  Sf  (Entropy  Production). 

in  the  following  paragraphs,  the  blade  force  and  entropy  production  terms 
are  discussed  while  discussion  of  the  mean  pressure  is  reserved  for  Section 
4.2.1. 

The  blade  force  Eg  of  Equation  2 represents  the  blade  force  acting  upon 
the  fluid.  The  blade  force  can  be  determined  through  reference  to  the  follow- 
ing sketch: 
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and  is  derived  from  the  Euler  Turbine  Equation.  There  is  an  unsteady  flow 
contribution  to  this  equation  which  is  discussed  in  Section  4.2.2.  The 
direction  of  the  lift  vector  is  assumed  to  be 
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which  is  based  upon  an  analogy  with  the  drag  coefficient  for  duct  flows.  It 
is  noted  that  the  drag  force  acts  in  the  direction  of  the  inlet  flow  Bi-  The 
term  Fp^  in  Equation  5 is  then  obtained  from  the  relation 
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where  the  numerator  transforms  the  force  (Fp)  in  the  flow  direction  to  a 
force  along  the  mean  chord  line  perpendicular  to  the  lift  vector  and  the 
denominator  transforms  that  resultant  force  to  an  axial  force. 

It  should  be  noted  that  in  steady  flow  a momentum  balance,  in  general, 
will  not  give  the  same  total-pressure  rise  per  stage  as  does  an  energy 
balance.  The  reason  for  this  difference  is  that  the  assumed  direction  of  the 
blade  lift  vector  is  not  exactly  the  arithmetic  average  (Equation  7)  of  the 
flow  angles.  Comparison  of  the  steady-state  momentum  and  energy  balance 
solutions  permits  the  determination  of  a small  "correction  angle"  which  can 
then  be  added  to  B<*>  to  give  the  proper  lift  direction.  Hence,  in  actuality 
the  lift  direction  angle  is  written  as 


Pc„  = 1/2  (B.  + B-)  + 8 (10) 
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Equation  7,  then,  is  replaced  by  Equation  10. 

The  term  Sp  in  Equation  3 represents  the  total  rate  of  irreversible 
conversion  of  mechanical  to  internal  energy  and,  in  the  case  of  this  model, 
represents  the  entropy  production  due  to  blade  row  losses.  It  can  be 
obtained  from  the  expression: 
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where  the  ideal  relative  total-pressure  ratio  which  accounts  for  the  change 
in  pitch  line  radius  from  the  entrance  of  a rotor  blade  row  to  its  exit  is 
written  as 


Mj  is  equal  to  the  ratio  of  the  blade- row- exit  pitch-line  wheel  speed  to  the 
inlet  relative  stagnation  velocity  of  sound  (ZirNr^aq^)  • In  the  case  of  a 
stator,  the  ideal  relative  total-pressure  ratio  is  equal  to  one.  The  actual 
relative  total-pressure  ratio  requires  knowledge  of  the  relative  total-pressure 
loss  coefficient  which  is  defined  as 
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Equation  13  can  be  rewritten  in  the  form 
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Hence,  Equations  12  and  14  when  substituted  into  Equation  11  provide  complete 
definition  of  the  entropy  production  term. 


As  might  be  expected,  the  input  to  the  program,  in  addition  to  physical  speed. 
Inlet  conditions  and  compressor  geometry,  requires  the  relative  total-pressure 
loss  coefficient  (to  ) for  each  blade  row,  the  tangent  of  the  deviation  angle 
(6)  for  the  rotors,  and  the  tangent  of  the  correction  angle  (Ec)  for  both  the 
rotors  and  stators.  The  relative  total-pressure  loss  coefficient  and  the 
deviation  angle  are  derived  from  stage-stacking  results  and  are  based  upon 
clean- inlet-flow  test  data  obtained  by  throttling  at  constant  speed.  The 
correction  angle  is  obtained  by  comparing  steady-flow  force  and  energy  balance 
solutions  on  a blade-row  basis.  These  parameters  can  be  represented  as  func- 
tions of  incidence  angle  and  are  input  to  the  program  in  this  manner.  These 
parameters,  in  conjunction  with  the  velocity  triangles  and  other  ancillary 
relations,  permit  the  determination  of  the  thermodynamics  of  the  fluid  at 
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■ nch  station.  Blade-free  volumes  are  treated  as  lossless  volumes  with  no 

ed  blade  force;  hence,  the  Fg  and  Sp  terms  of  Equations  2 and  3 are 

: ien finally  zero. 

line  dependent  solution  of  the  system  of  equations  (Equations  1 through 
3 md  the  relations  for  Eg,  P^,  and  Sp)  that  comprise  the  dynamic  digital 
.■■.impression  component  model  is  effected  through  a Taylor  series  which 

■ stablishes  the  values  of  the  three  independent  volume-averaged  variables  at 
the  next  increment  in  time.  In  the  case  of  this  model  and  with  references  to 
the  left  hand  side  of  Equations  1,  2,  and  3,  the  variables  p,  W,  and  ps  are 
ones  for  which  a solution  is  sought.  Solution  is  now  straight-forward  and 
will  be  illustrated  for  one  variable  - the  volume-averaged  density.  Consider- 
ing that  this  method  is  applicable  to  any  volume,  the  subscript  "k",  indicating 

ii  volume  will  be  dropped.  The  Taylor  series  for  volume-averaged 
1 rr.-  t to  second  order  can  be  written  as: 


p ( t+  ’.t) 


P (t) 


+ 9P  (O 
3t 


At  + 


(t)  (At)^ 

2 2 


(15) 


where : 

p(t)  is  established  by  the  initial  conditions  or  from  the  previous  time  step. 


3p(t)  = 1 
3t  V 


(W. 


W.+1)  from  Equation  1 and  differentiating 


Equation  1 with  respect  to  time  yields: 


(t) 

2 


(16) 


Examination  of  Equation  16  reveals  that  the  right-hand  side  is  composed 
of  derivatives  of  station  values  of  flow  with  respect  to  time.  Since  Equation 
2 will  supply  only  the  derivatives  of  the  volume-averaged  flow  with  respect 
to  time,  use  of  an  interpolation  scheme  for  obtaining  station  values  from 
volume-averaged  values  will  permit  Equation  16  to  be  solved  for  the  second 
partial  derivative  of  volume-averaged  density  with  respect  to  time.  Equa- 
tions 1 and  16  then  can  be  substituted  into  Equation  15  to  obtain  the  estimate 
of  the  volume  averaged  density  correct  to  second  order  at  the  next  increment 
in  time.  Equation  16  implies  that  first  derivatives  with  respect  to  time  of 
a large  number  of  terms  (e.g.,  Fg,  P^,  and  Sp)  will  be  required.  Although 
these  expansions  are  lengthy,  they  can  be  derived  in  a straightforward  manner 
and  will  not  be  reproduced  here.  Similarly,  this  technique  can  be  used  for 
the  remaining  two  variables  (W  and  ps)  and  can  be  continued  from  one  time 
step  to  the  next  for  the  desired  number  of  time  steps. 
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This  calculational  scheme  is  numerically  stable,  and  (to  date)  anomalous 
behavior  has  occurred  only  when  a physical  aerodynamic  instability  in  the  flow 
would  be  expected  to  occur.  There  are  two  related  aspects  of  the  errors  intro- 
duced by  truncating  the  Taylor  series:  1)  error  in  estimation  of  a parameter 
at  a given  time  and  2)  buildup  of  error  with  time.  The  error  in  a parameter 
at  a given  time  introduced  by  truncation  is  given  in  Lagrangian  form  (Refer- 
ence 11)  by 


R2  = 


(At)3 

3! 


. 3 , 

3 p (t) 

"at  3 


t<  t < t+At 


(17) 


This  error  is  negligibly  small  for  typical  time  increments  on  the  order  of 
1 0— ^ to  10“ ^ seconds.  The  error  introduced  by  truncating  the  Taylor  series 
after  second  order  does  not  appear  to  give  rise  to  an  error  buildup  in  the. 
calculated  parameters  as  time  increases.  Rather  these  errors  result  in  con- 
tributions to  the  time  derivatives  in  the  macrobalances  which  then  act  as 
restoring  forces.  This  behavior  has  been  noted  during  simulation  of  steady- 
<tate  operating  conditions. 

Now  the  calculational  technique  utilized  in  the  Dynamic  Digital  Blade 
Row  Compression  Component  Stability  Model  can  be  discussed  and  is  illus- 
trated in  Figure  10  in  block  diagram  format.  The  use  of  the  blade  row  build- 
ing block  concept  allowed  construction  of  a generalized  model  which  is  indepen- 
dent of  the  particular  compression  component  being  simulated. 

Block  I is  a statement  of  the  required  dependent  variable  information, 
that  is,  volume-averaged  density,  flow,  and  entropy.  These  variables  are 
available  from  either  a steady-state  (SS)  initialization  or  a previous  time 
step  of  a time-dependent  (TD)  analysis. 

Block  II  presents  the  macrobalances  in  the  form  they  are  used  in  the 
analysis.  The  variables  on  the  right-hand  side  of  the  equations  are 
station-value  properties.  Knowledge  of  these  parameters  allows  the  first 
time  derivatives  of  the  volume-averaged  properties  to  be  calculated.  How- 
ever, as  stated  in  Block  I,  only  volume-averaged  quantities  are  available 
at  the  beginning  of  each  time  step.  Therefore,  Block  II  illustrates  that  it 
is  necessary  to  interpolate  between  volume-averaged  parameters  in  order  to 
obtain  station-value  properties. 


In  the  case  of  blade-free  volumes,  where  no  blade  forces  act  or  entropy 
production  takes  place,  it  is  only  necessary  to  calculate  station  axial 
velocity  in  order  to  evaluate  the  macrobalances.  As  shown  in  the  lower 
branch  of  Block  IV,  the  assumption  of  constant  absolute  flow  angle  across 
the  volume  is  made.  Total  pressure,  total  temperature  and  other  desirable 
parameters  are  also  calculated  at  this  point.  A special  case  of  the  blade 
free-volume  calculations  is  the  imposition  of  the  boundary  conditions.  At 
the  model  inlet,  constant  total  pressure  and  total  temperature  as  well  as 
constant  entropy  are  maintained.  At  the  exit  of  a model  of  a turbojet  com- 
pressor and  burner,  a specified  exit-flow-function  boundary  condition  is 
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imposed.  This  boundary  condition  was  derived  from  the  assumption  that  a 
fictitious,  zero-length  choked  nozzle  existed  at  the  burner  exit.  The 
fictitious,  zero-length  choked  exit  condition  implies  the  upstream  flow 
function  at  the  burner  exit  is  specified  and  can  be  either  held  constant  or 
changed  at  some  rate  to  simulate  a constant-speed  throttling  process. 

As  indicated  in  the  upper  branch  of  Block  IV,  the  presence  of  a bladed 
volume  requires  the  net  axial  blade  force,  entropy  production,  and  station 
axial  velocities  to  be  calculated.  Calculation  of  the  net  axial  force  and 
entropy  production  terms  require  knowledge  of  the  loss  coefficient  and  devi- 
ation angle.  This  information  is  available  as  polynomial  representations 
which  are  functions  of  incidence  angle.  Stationary  blade  rows  may  be  assumed 
to  be  lossless  and  are  assumed  to  have  constant  deviation  angles. 

Once  the  flow  conditions  at  the  stations  are  completely  described, 
various  quantities  of  interest  can  be  calculated  such  as  stage  coefficients, 
diffusion  factors,  etc.  With  all  the  necessary  quantities  on  the  right  hand 
side  of  the  macrobalances  available,  the  first  time  derivatives  of  the  volume- 
averaged  properties  can  be  calculated  as  indicated  in  Block  V using  the  macro- 
balances of  Block  II. 

Expressions  for  the  second  time  derivatives  of  the  volume  averaged 
quantities  are  obtained  from  differentiating  the  macrobalances  with  respect 
to  time.  Analytical  expressions  for  the  time  derivatives  of  the  station 
properties  can  be  evaluated  by  interpolating  between  volumes  and  through 
use  of  the  macrobalances.  Thus,  as  indicated  in  Block  VI,  the  second  time 
derivatives  of  volume-averaged  properties  can  be  calculated. 

This  procedure  for  calculating  station  properties  and  evaluating  first 
and  second  time  derivatives  of  the  volume-averaged  properties  can  be  carried 
out  for  any  number  and  types  of  volumes  (Block  VII)  and  is  not  dependent  on 
the  particular  geometry  being  modeled.  Once  these  calculations  are  carried 
out  for  all  the  volumes,  the  solution  can  be  advanced  to  the  next  time  step 
through  use  of  the  second-order  Taylor's  series  approximations,  Block  VIII. 

As  specified  in  Block  IX,  the  technique  can  be  repeated  for  as  many  time 
steps  as  required  by  the  event  being  simulated. 


A . 2 DIGITAL  COMPRESSION  COMPONENT  MODEL  REFINEMENTS 

During  this  contractual  effort,  a number  of  refinements  were  introduced 
to  the  basic  model  as  reviewed  in  the  previous  sections.  These  refinements 
were  introduced  to  improve  the  generality  of  the  model,  its  accuracy,  and  to 
handle  inherently  unsteady  flows  correctly,  since  prior  to  this  study,  the 
model  had  been  exercised  in  detail  only  for  slow  throttling  transients  ot  a 
compressor.  The  following  paragraphs  detail  a generalized  method  for  handling 
the  free-volume  mean-pressure  term.  In  addition,  the  formulation  for  the 
unsteady  contribution  to  the  tangential  component  of  the  blade  force  is  given 
as  well  as  a method  for  handling  the  lift  force,  loss  coefficient,  deviation 
angle,  and  lift  direction  correction  angle  in  a way  such  that  at  high 
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frequencies  it  is  not  necessary  to  assume  that  these  parameters  adjust  instan- 
taneously to  a new  value  of  incidence  angle. 


4.2.1  Free-Volume  Mean  Pressure 

The  mean  pressure  Pm  of  Equation  2 represents  the  integral  of  the  pres- 
sure distribution  acting  over  the  lateral  surface  area  of  a volume.  Its 
relationship  to  the  outer  terms  in  the  momentum  equation  is  shown  in  the 
sketch  below 


where : 


(18) 


The  expression  that  was  used  for  the  free  volume  mean  pressure  in  the 
Reference  4 study  was  based  upon  an  empirical  relationship  involving  a scale 
factor.  For  high  Mach  numbers  or  for  large  area  changes  across  a volume, 
this  expression  was  too  inaccurate  and  led  to  numerical  stability  problems. 
Therefore,  an  effort  was  undertaken  to  improve  the  accuracy  of  Pm,  its 
generality,  and  to  provide  it  with  some  analytical  basis.  One  aspect  that 
must  be  considered  in  developing  a relation  for  PM  is  that  it  must  be 
independent  of  the  momentum  equation  (Reference  2)  in  order  that  redundancy 
is  not  encountered. 

As  a starting  place  for  our  derivation,  we  begin  with  Equation  18 
written  in  the  form 
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(19) 


which  can  be  done  if  the  total  pressure  and  total  temperature  across  a free 
volume  are  constant.  During  steady  flow,  the  total  pressure  and  total  tem- 
perature across  a lossless  free  volume  will  remain  constant  and  hence,  the 
flow  will  be  isentropic.  Therefore,  the  isentropic  relations  for  P/P-p  and 
A/A  can  be  written  as 

P/PT  = (1  + M2)  ~ y-1  (20) 


_Y±L_ 

2(y-1) 

A/A*  = M 

The  differentiation  and  integration  tasks  implied  by  Equation  19  can  be  sig- 
nificantly reduced  by  curve  fitting  Equations  20  and  21  with  polynomials  as 
follows : 

P/P,r  = Aj  + B M + C^M2  + D^3  + F^M4  (22) 

A/A*  = A2  + B2M  + C2M2  + D2M3  + E2M4  (23) 

where  the  values  of  the  coefficients  are  given  in  Table  2. 

Table  2.  Isentropic  Relationships  Polynomial  Coefficients 


Coef  f icient 

P/PT 

A/A* 

A 

0.9999810 

1.4648536 

B 

0.0018162 

-1.0785448 

C 

-0.7229570 

0.3238214 

D 

0.0967896 

-0.0407312 

F. 

0. 172 J 226 

0.0017179 

2(1  + M2) 

~ y+t 


(21) 
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Equation  23  can  be  differentiated  and  substituted  into  Equation  19  along 
with  Equation  22.  The  integrations  in  Mach  number  can  easily  be  carried 
out  to  give 

A (M  -M  ) + B (M^V)/2  + C,(M3-M3)/3  + D_(M4V)/4  + E.(M^-mJ)/5 

P /p  = — — i i— t — i 2 — £ — ± i — £ k ill  c24 ■) 

M T 22  33  44 

B2(M2-M1)  + C2(M2-M1)  + D2(M2-Mp  + E2(M2-M1) 


where 


A3  " A1B2 


B3  = BjB.,  + 2A]C2 


C3  = C1B2  + 2B1C2  + 3A1°2 


D3  = DlB2  + 2ClC2  + 3BlD2  + 4AlE2 


E3  = ElB?  + 2DlC2  + 3ClD2  + 4BlE2 


(25) 

(26) 

(27) 

(28) 
(29) 


During  the  reported  effort,  an  approximation  to  Equation  24  was  used 
which  can  be  derived  from  the  relation 


P /P 
M T 


r2 

J P/PT  dM 

r 

j J dM 


(30) 


Upon  carrying  out  the  integration  of  Equation  30,  the  following  expression 
is  obtained 


P /P  = 
M'  T 


A1(M2-M1)  + Bp  /2  + Ci(M2"M1)/3  + D1(M2_M1)/4 

R__ 


(31) 


which  is  seen  to  be  the  lead  terms  of  Equation  24.  Equation  31  matches 
Equation  24  within  a few  percent  in  the  range  of  interest. 


Although  the  relations  given  by  Equations  24  and  31  were  derived  for 
steady  flow,  they  are  assumed  to  apply  in  unsteady  flows  where  the  wavelength 
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of  the  disturbance  is  considerably  greater  than  the  length  of  a volume. 
Further,  since  in  unsteady  flows  the  total  pressure  is  not  generally  con- 
stant, it  is  necessary  to  assume  that  the  total  pressure  P-p  of  Equations  24 
and  31  is  equal  to  the  average  of  the  inlet-  and  exit-total  pressures. 


4.2.2  Tangential  Blade  Force  in  Unsteady  Flows 

The  tangential  blade  force  is  determined  from  a generalized  expression 
for  the  conservation  of  angular  momentum.  As  our  starting  point,  we  begin 
with  the  moment  of  linear  momentum  equation  (Reference  12)  f or  a stationary 
control  volume  in  the  form 


l (F  x r) 


(C 


x r)  C • djp  + ~ j'p  CxrdR  j 

D 


(32) 


The  first  term  of  the  right-hand  side  of  Equation  32  represents  the  net  change 
in  flux  of  angular  momentum  across  the  control  surface  S and  the  second  term 
represents  the  rate  of  change  of  angular  momentum  of  the  fluid  within  the 
region  R.  Only  the  angular  momentum  components  which  contribute  to  the  net 
torque  are  considered.  Hence,  only  tangential  velocity  components  will  con- 
tribute, since  the  radial  velocity  acts  through  the  origin. 


In  steady  flow.  Equation  32  reduces  to  the  well  known  Euler  turbine 
relationship  in  the  following  manner.  The  fluid  is  assumed  to  be  homogenous 
within  the  region  and  uniform  at  the  boundaries  of  the  surface  at  two  geo- 
metrically defined  locations.  Then,  the  factor  pV  • d5p  is  simply  equal  to 
the  differential  of  flow  rate  dW.  Hence,  the  integral 


y*p  (C  x r)  C • dS  = / (C  x r 


) dW 


Considering  the  tangential  component  of  velocity,  then 


p (C  x r)  C • dS  = r2  Cu  W2  - r^ 


and  for  steady  flow  Equation  32  becomes 


FT  r ■ i„  <r2  Cu2  “2 


or 


r.  C W ) 

1 Uj  1 


FT  = <r2  Cu2  W2 


ri  Cu1  "A  (ri  + r2)  / 2j 


(33) 


(34) 


(35) 


(36) 
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Now  let  us  consider  the  second  term  on  the  right-hand  side  of  Equation 
32.  The  storage  of  tangential  angular  momentum  can  be  written  as 


dt  /p  ^ x £ dR  = ^ f P cu 


r dR  = (pC  r V) 
d t u 


/ P C x r dR  = Vr 


TT  (p  c ) 
dt  u 


Therefore,  Equation  32  can  now  be  written  in  the  form 


F_  r = ~ (r9  C W - r.  C W, ) + Vr 

f gQ  2 u2  2 1 u1  1 dt  i 


2 (r  C W - r C W ) 

i.  1 , ui  .i  „ v 4^) 

go  rj  ♦ r2  dt  k 


Note  that 


(C  + C ) 

d (pCu)  k „ d_  _^1 ^2 

dt  dt  P 2 


Pk  U1  U2 

— — ( + -) 

2 ( dt  dt  ' 


+ C„  + C„ 

U1  U2  ^k 


Inspection  of  Equation  41  with  respect  to  the  manner  in  which  compu- 
tations are  carried  out  in  the  model  shows  that  all  terms  on  the  right-hand 
side  are  known  with  the  exception  of  dC^/dt.  Solution  of  the  momentum 
equation  given  by  Equation  2 for  the  time  rate  of  change  of  the  volume 
averaged  flow  in  the  k-th  volume  requires  knowledge  of  the  derivative  with 
respect  to  time  of  the  absolute  tangential  velocity  at  the  exit  of  the 
volume.  But,  the  time  rate  of  change  of  the  absolute  tangential  velocity 
at  the  exit  cannot  be  determined  unless  the  time  rate  of  change  of  the  volume 
averaged  flow  has  been  determined.  This  impasse  which  occurs  within  a time 
step  can  be  avoided  by  using  information  from  the  previous  time  step  in  the 
following  manner. 
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Since  CU2  (t  - At)  and  (^(t)  are  known  values  prior  to  making  the  cal- 
culation for  the  time  rate  of  change  of  the  absolute  tangential  velocity,  the 
desired  derivative  can  be  approximated  as  follows: 


d C (t) 

ll2 

dt 


and  is  updated  at  each  time  step.  Detailed  calculations  show  that  compu- 
tation of  the  derivative  in  this  manner  provides  an  accurate  approximation 
of  the  precise  value  of  dC^/dt. 
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4.2.3  Unsteady  Blade  Characteristics 

The  lift  coefficient  (Reference  13)  for  thin,  non-fluttering,  non- 
translating isolated  airfoils  immersed  in  an  oscillating  flow  can  be 
written  as 


— \ = *C(k)a  + | 11  + C(k)  ] a (43) 

pU  c/2  » 

CO 


where  the  term  proportional  to  a gives  the  contribution  of  lift  proportional 
to  the  mean  angle  of  attack  and  the  term  proportional  to  a gives  the  contri- 
bution of  lift  proportional  to  the  change  in  circulation  about  the  airfoil 
with  time. 

An  equivalent  expression  for  a cascade  does  not  exist.  But  an  approxi- 
mation can  be  formulated  which  treats  the  essence  of  the  problem.  The  work 
of  Melick  (Reference  14)  suggested  an  approach  which  established  a relation- 
ship between  the  instantaneous  angle  of  attack  of  a blade  in  a cascade  and  an 
effective  angle  of  attack  of  a blade  which  is  associated  with  the  actual  values 
of  the  blade  lift  coefficients.  This  approach  has  subseqeuntly  been  drawn 
upon  by  Kimzey  (Reference  15)  for  treating  the  effect  of  varying  angular 
extent  1/rev  circumferential  distortions. 
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It  is  postulated  that  the  effective  angle  of  attack  of  an  airfoil  in 
unsteady  flow  can  be  related  to  the  instantaneous  angle  of  attack  of  the  air- 
foil  through  a differential  equation  of  the  following  form: 


d(a  -a  ) 

--4-2-  + 

dt 


( tx  u ) 
E o 


(ci  -a  ) 
1 o 


(46) 


where 


aE  = 

effective  angl 

e of  attack 

«i  = 

instantaneous 

angle  of 

attack 

(X  = 

o 

mean  angle  of 

attack 

T = 

time  constant 

for  n-th 

blade 

It  i noted  t dal  in  the  limit  ol  steady  flow,  the  effective  angle  of  attack 
. equal  to  the  instantaneous  angle  of  attack.  Given  that  the  instantaneous 
angle  of  attack  is  known,  it  is  possible  to  find  the  effective  angle  of 
ittm  k if  the  time  constant  for  each  blade  row  is  known. 


The  work  of  Schorr  and  Reddy  (Reference  8)  provides  predictions  for  a 
cascade  with  oscillating,  inviscid,  incompressible  inlet  flow  from  which  blade 
row  time  constants  can  be  derived.  Their  results  in  terms  of  a nondimensional 
dynamic  lift  coefficient  as  a function  of  reduced  frequency  are  reproduced 
with  some  slight  changes  in  nomenclature  in  Figure  11.  Further,  their  work 
shows  that  the  dynamic  lift  coefficient  is  a stronger  function  of  solidity 
than  stagger  angle  in  the  region  of  interest  (oi*  _<  3.0).  Equation  46  can 
be  solved  for  the  Schorr  and  Reddy  parameter  by  assuming  aj  - a0  = Aa0  sin  uot. 
It  then  becomes  necessary  to  choose  the  value  of  the  time  constant  Tn  which 
causes  the  solution  of  Equation  47  and  the  results  of  Schorr  and  Reddy 
to  match. 


d (a  -a  ) 
F.  o 


(a  -a  ) = A a sin  wt 
E o o 


(47) 


1 i particular  solution  to  this  equation  can  be  written  as 


Au 


a -a 
E o 


sin  (cot  - 0) 


(48) 


V l+(t u) 


and  the  absolute  magnitude  of  the  dynamic  fluctuations  can  be  written  as 


3<> 


I 


(49) 


Aa 

o 


\/l+(i  w)2 
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If  the  angle  of  attack  Is  assumed  proportional  to  the  lift  coefficient  then 
we  obtain  the  ordinate  of  Schorr  and  Reddy 


1C. 


L | AMP 
L SS 


where 


T W 

n 


and 


1 


v l+(x  w)2 
n 


c c 

-jj-  2 7T  f = C-W* 
R 


(50) 


(51) 


C = Blade  Chord 

U„  = Relative  Velocity 
K 

f = Frequency 

w*  = Reduced  Frequency 

= Constant  Dependent  on  Blade  Row  Solidity 

Figure  11  has  been  replotted  in  Figure  12  and  extrapolated  curves  for 
solidities  greater  than  1.25  have  been  added.  These  curves  have  been  curve 
fit  using  Equations  50  and  51  to  determine  the  constant  C3.  This  constant  is 
shown  in  Figure  13  as  a function  of  solidity. 

The  transformation  of  these  results  for  use  in  the  Dynamic  Digital  Blade 
Row  Compression  Component  Stability  Model  is  straight  forward.  Equation  46 
can  be  written  in  terms  of  incidence  angle  by  noting  for  circular  arc  air- 
foils 


a = 1 + 0/2 


(52) 


where  0 is  the  blade  camber  angle  as 


T 

n 


d(iE-V 


dt 


+ (iE-V 


(irV 


(53) 
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Effective  Incidence  Angle 
Mean  Angle  of  Incidence  Angle 
Instantaneous  Incidence  Angle 

C3 


Since  the  model  computes  the  instantaneous  incidence  angle  for  each  blade  row, 
then  given  the  solidity,  blade  chord,  and  relative  velocity  for  a blade  row, 
the  blade-row  time  constant  is  defined  and  Equation  53  can  be  solved  for  the 
effective  incidence  angle  through  use  of  a second-order  Tay lor-ser ies  expan- 
sion. In  effect,  the  relative  total-pressure  loss  coefficients,  deviation 
angles,  lift  direction,  and  lift  direction  correlation  angles  become  func- 
tions of  the  effective  incidence  angle  rather  than  the  instantaneous  incidence 
angle  as  was  assumed  in  the  original  model  formulation. 

The  assumption  inherent  to  this  approach  is  contained  in  assuming  that 
Schorr  and  Reddy's  work  for  incompressible  flow  is  also  applicable  to  com- 
pressible flow.  In  actuality,  it  is  anticipated  that  detailed  comparisons 
between  model  predictions  and  test  data  will  result  in  the  modification  of 
the  relationship  between  solidity  and  the  constant  C3  of  Figure  13.  The 
effect  of  including  the  aspects  of  unsteady  blade  characteristics  as  embodied 
in  Equation  53  is  discussed  in  Section  5.5. 


i) . 3 DIGITAL  MODEL  OF  PLANAR  PRESSURE  PULSE  GENERATOR 

The  planar  pressure  pulse  generator  is  a siren-like  device  with  a stator 
and  rotor  positioned  in  such  a manner  that  the  mean  effective  area  is  modu- 
lated as  the  rotor  rotates  producing  variations  in  the  total-pressure  loss  of 
the  air  drawn  through  the  rotor-stator  combination  by  a downstream  energy 
source.  In  the  case  of  the  present  study,  the  energy  is  supplied  by  a two- 
stage  fan.  The  P^G  area  is  essentially  choked  at  all  times  and  the  total- 
pressure  loss  results  primarily  from  shock  losses  as  the  air  upstream  of  the 
P^G  attempts  to  adjust  to  the  downstream  Mach  number  boundary  condition. 

Some  viscous  losses  are  present,  but  for  the  flow  rates  and  planar  wave 
amplitudes  of  interest,  the  shock  losses  predominate.  The  amplitude  of  the 
pressure  pulses  are  physically  controlled  by  modulating  the  separation  between 
the  rotor  and  stator  since  the  stator  translates  axially.  This  changes  the 
effective  area  (both  in  the  mean  and  at  the  extrema)  and,  therefore,  varies 
the  amplitude  of  the  total-pressure  loss. 

The  P^G  simulation  approximates  the  aerodynamics  of  the  P^G  operation  as 
follows:  A steady-state  condition  satisfying  the  "mean"  conditions  at  the 

P^G  during  operation  is  established  using  the  mean  steady-state  recovery 
across  the  P^G  as  measured  during  the  test.  This  recovery  implies  an  effective 
choked  area  at  the  P^G  which  can  be  computed.  The  mean  area  becomes  the  base 
for  the  amplitude  of  the  area  fluctuation.  The  amplitude  is  controlled  by  the 
parameter  "amp"  as  shown  in  the  following  equation. 


[1  + amp  • sin (2nNt/60) ] 


choked 


The  choked  area  can  be  related  to  the  upstream  duct  A*/A(juct  and  provides 
a fundamental  boundary  condition  to  the  P^C  problem.  As  much  of  the  existing 
General  Electric  Dynamic  Digital  Blade  Row  Compression  Component  Stability 
Model  logic  was  used  as  possible  and  the  P-^G  simulation  was  constructed  pri- 
marily as  a subroutine  to  the  main  program.  With  reference  to  the  model  as 
described  in  Section  4.1  and  schematically  described  in  the  block  diagram  of 
Figure  10,  the  P^G  subroutine  was  treated  as  a block  parallel  to  the  free 
volume  and  blade  row  subroutines  in  Step  IV.  All  necessary  derivatives  are 
computed  as  additions  to  the  blocks  of  Steps  V and  VI. 


4.4  DIGITAL  MODEL  OF  TEST  INSTALLATION 

The  system  geometry,  the  two-stage  fan-blade-row  characteristics,  and 
the  boundary  conditions  imposed  on  the  model  are  discussed  in  the  following 
paragraphs . 


4.4.1  System  Geometrv 


The  geometry  of  the  test  configuration  was  shown  previously  in  Figure  2. 
The  entire  test  set-up  from  the  bellmouth  entrance  to  the  bypass  discharge 
valve  exit  was  simulated  by  the  system  model.  Table  3 defines  the  individual 
characteristics  of  the  69  volumes  used  to  model  the  overall  facility  and  the 
area  distribution  of  the  test  setup  as  used  in  the  model  is  shown  in  Figure 
14.  Each  blade  row  of  the  fan  is  represented  by  one  volume  such  that  the 
minimum  volume  length  is  2.54  inches  and  is  dictated  by  the  axial  length  of 
the  second  stator.  The  maximum  and  most  common  volume  length  is  4 inches  and 
is  used  to  represent  the  blade-free  volumes  in  the  system.  The  IGV  was 
represented  by  two  distinct  volumes;  one  for  the  forward-fixed  portion  of  the 
vane  and  one  for  the  movable  trailing  edge.  This  representation  of  the  IGV 
was  accomplished  without  compromising  the  blade  row  aerodynamics  and  per- 
mitted maintaining  the  length  of  all  volumes  equal  to  or  less  than  4 inches. 
If  the  IGV  were  not  represented  by  two  volumes,  its  length  would  have  been 
7.7  inches  and  the  model  frequency  response  would  have  been  sacrificed.  The 
total  length  of  the  modeled  system  was  266  inches. 

The  physical  characteristics  of  the  fan  geometry,  the  metal  angles, 
pitch-line  radii,  blockages,  areas,  and  solidities,  are  also  defined  for  each 
blade  row  in  Table  3.  In  addition,  the  exit  air  angle  for  each  stator  is 
given.  Table  4 shows  the  IGV  exit  metal  and  air  angle  schedules  as  a func- 
tion of  corrected  speed  (100%  N//0  = 7556  rpm) . 


Figure  lk.  System  Area  Distribution 
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The  location  of  instrumentation  is  shown  also  in  Figure  2 and  Table  3. 
Instrumentation  locations  that  were  used  during  the  test  were  duplicated  in 
the  model  setup  to  provide  a direct  comparison  between  model  results  and  test 
data.  The  parameters  chosen  for  this  comparison  are  defined  in  Table  5 and 
were  consistent  with  the  objectives  of  verifying  model  fidelity  as  well  as 
documenting  the  fan  aerodynamic  response. 

Table 

5. 

Model  Parameters  Used  in  Validation  Study. 

Parameter 

Description 

P3GPS 

Static  pressure  between  the  forward  boundary 
(bellmouth  inlet)  and  the  P^G  plane. 

AXPS2  1 
AXPS3  f 

Static  pressures  in  the  duct  between  the  P^G 
and  the  fan  inlet. 

PTIGV 

Total  pressure  at  the  IGV  plane. 

WCA 

Approximate  calculated  airflow  at  the  IGV  plane. 

PS14 

Static  pressure  at  the  fan  exit. 

FDPSC 

Static  pressure  in  the  fan  duct  just  down-stream 
of  the  splitter. 

FDPSD 

Static  pressure  close  to  the  bypass  discharge 
in  the  fan  discharge  duct. 
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The  tan  duct  discharge*  area  was  divided  at  the  core  entrance  plane  in 
the  physical  test  set-up.  For  the  purposes  of  the  model,  the  two  flow  areas 
(core  dut  and  bypass  duct)  were  summed  together  and  represented  b/  a single 
t low  area. 


4.4.2  Iwo-Stage  Fan  Characteristics 

The  performance  characteristics  of  the  two-stage  fan  were  established 
during  steady-state  component  testing.  During  this  type  of  testing,  inter- 
stage aerodynamic  information  was  obtained  from  which  the  performance  of  the 
ndividual  stages  was  established  and  the  performance  of  the  blade  rows  was 
inferred.  Internal  fan  steady-state  pressures  and  temperatures  were  measured 
is  lb  fin  was  throttled  from  its  lowest  operating  point  until  surge  occurred. 

;.d  up  n standard  compressor  design  techniques,  this  data  was  represented 
in  two  different  forms.  The  first  was  a non-dimensional  representation  of 
the  stage  characteristics,  showing  the  pressure  coefficient  (0'),  the  work 
left;,  ient  ('p)  and  stage  efficiency  fit)  as  functions  of  flow  coefficient 
t)  for  each  stage.  The  positive  feature  of  this  format  was  that  the  form  is 
amt nable  to  generalization  and  that  the  behavior  of  a stage  may  be  assessed 
wif'i  reference  to  stage  characteristics  established  for  other  compression 
omponents.  This  form  Is  useful  for  stage  performance  representation,  but 
not  for  blade-row  by  blade-row  characteristics  which  are  the  basic  input  to 
the  model  simulation.  Blade-row  performance  characteristics  take  the  form  of 
blade-row  relative  total-pressure  loss  coefficients  (w')  and  deviation  angles 
(6)  as  a function  of  air  incidence  angle  (i).  In  this  program,  all  losses 
were  issigned  to  the  rotor  and  it  was  assumed  that  only  lossless  turning  of 
the  flow  occurs  in  the  stators.  This  assumption  is  necessary  if  high-quality 
'est  data  are  not  obtained  at  each  blade  row  or  design  stator  loss  coefficients 
are  not  available.  Although,  only  the  blade  characteristics  form  of  per- 
formance representation  is  directly  applicable  to  the  analytical  model, 
the  stage  characteristics  are  useful  in  analyzing  the  credibility  of  the 
assumptions  that  go  into  the  blade-row  by  blade-row  results. 

The  rotor  loss  coefficients  and  the  deviation  angles  as  functions  of 
incidence  angle  (tan  i)  are  tabulated  in  Appendix  B for  both  rotors  at  the 
•wo  test  speeds  of  interest,  80  and  100%,  respectively.  These  character- 
istics represent  blade-row  performance  up  to  steady-state  surge.  However, 
due  to  the  fundamental  nature  of  the  disturbances  that  were  to  be  imposed  on 
the  model,  it  was  known  that  the  incidence  angles  encountered  during  model 
operation  would  instantaneously  exceed  the  steady-state  surge  limit.  There- 
fore, an  analysis  was  conducted  of  the  extrapolation  of  the  steady-state  data 
into  the  post  steady-state  surge  regime.  Numerous  extrapolations  were  attempted 
and  analyzed  using  the  non-dimensional  - ip , \p  ’ , q format  to  aid  in  the 
interpretation  of  results  in  this  uncharted  region.  At  the  conclusion  of 
this  study,  it  was  determined  that  the  most  simple  of  data  extrapolations 
worked  best,  a straight  line  tangent  to  the  data  at  the  end  points.  Any 
other  technique  produced  an  unrealistic  "folding  back"  of  the  characteristics 
that  did  not  appear  to  be  consistent  with  the  aerodynamics  of  the  compression 
component.  Subsequent  model  results  have  shown  the  importance  of  the 


i 


extrapolated  characteristics,  but  the  detailed  shape  in  the  steady-state 
post-surge  region  still  cannot  be  established  with  the  existing  model  and 
existing  test  results. 

The  blade  row  performance  characteristics  are  represented  by  a fourth 
order  polynomial  curve  fit  of  the  steady-state  data  with  linear  extrapola- 
tions at  the  end  points.  Figure  15  illustrates  this  technique  showing  the 
curve  through  the  data  and  the  end  points  of  the  polynomial  representation 
beyond  which  a linear  extraploation  of  the  curve  at  a slope  tangent  to  the 
end  point  is  used.  Plots  of  the  characteristics  and  tabulations  of  the  poly- 
nomial coefficients  as  used  in  the  model  are  contained  in  Appendix  B. 

It  should  be  noted  that  all  of  the  results  obtained  from  the  model 
during  this  study  were  obtained  using  blade  characteristics  that  were  ob- 
tained from  undeteriorated  performance  test  results.  However,  the  test 
results  against  which  the  model  results  are  being  compared,  inc  ude  pro- 
gressive amounts  of  deterioration  which  culminate  in  the  deteriorated  100% 
corrected  speed  line  as  was  shown  in  Figure  7. 

One  item  of  empirical  information  is  used  to  represent  the  test  data  in 
the  model,  the  lift  direction  correction  angle.  This  parameter  was  discussed 
in  Section  4.1  and  the  curves  of  this  parameter  and  tabulations  of  the  poly- 
nomial coefficients  are  given  in  Appendix  B. 


4.4.3  Boundary  Conditions 

The  system  model  is  bounded  by  two  distinctly  different  conditions  that 
are  intended  to  approximate  the  physics  of  the  environment  at  their  respec- 
tive locations. 

The  forward  boundary  position  is  geometrically  coincident  with  the  bell- 
mouth  inlet  plane.  Constant  entrance  total  pressure  and  total  temperature 
equal  to  the  measured  test  values  were  assumed  for  the  inlet  plane  aero- 
dynamic boundary  condition.  Subsequent  analysis  of  this  assumption  has  indi- 
cated that  the  bellmouth  entrance  conditions  are  not  exactly  "plenum-like," 
but,  are  sufficiently  close  for  all  practical  purposes.  Other  considerations 
which  encouraged  the  use  of  this  boundary  condition  were  the  sheer  size  of 
the  system  model  in  its  basic  state  which  necessarily  restricts  the  construc- 
tion of  a larger  plenum  forward  of  the  bellmouth  and  the  fact  that  the  P^G 
downstream  of  the  inlet  provides  a very  hard  boundary  condition  of  its  own, 
decoupling  the  fan  from  the  inlet  representation. 

The  exit  boundary  condition  was  established  at  the  bypass  discharge 
valve  location.  A constant  static  pressure  is  established  in  the  model  at 
this  point  and  is  felt  to  be  consistent  with  the  actual  test  conditions  since 
there  are  no  choke  points  or  even  high  Mach  numbers  downstream  of  the  fan 
exit . 
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4.5  METHOD  OF  MODEL  OPERATION 
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The  procedure  for  establishing  conditions  in  the  dynamic  model  similar 
to  test  conditions  evolved  as  the  simulations  were  attempted  until  the  model 
operational  technique  became  remarkably  similar  to  the  test  method.  As 
previously  described,  the  test  procedure  called  for  establishing  the  fan  at  a 
relatively  low  level  operating  point  on  the  chosen  speed  line  and  then  bring- 
ing the  P^G  to  the  appropriate  frequency  before  decreasing  the  rotor-to- 
stator  spacing  to  increase  the  planar  wave  amplitude.  Similarly,  the  model 
enters  the  computation  at  the  chosen  fan  speed  and  a near  surge  operating 
point.  The  initial  operating  point  was  not  always  commensurate  with  the  test 
reading  in  the  interests  of  conserving  computation  time  and  as  a result,  some 
dissimilarities  in  mean  point  behavior  occurred  between  model  and  test  results. 
However,  a mean  operating  point  was  established  on  the  model  consistent  with 
the  mean  total-pressure  loss  across  the  P^G  at  the  desired  rotor-to-stator 
(R-T-S)  spacing.  The  area  fluctuation  about  this  mean  condition  was  then 
increased  slowly  until  the  IGV  total-pressure  fluctuation  of  the  model  matched 
the  measured  value.  Behavior  of  the  model  P^G  representation  with  respect  to 
amplitude  versus  frequency  was  not  always  predictable  due  to  standing  wave 
patterns.  As  a result,  several  attempts  to  match  model  and  test  amplitudes 
were  required  for  some  cases.  The  procedure  of  ramping  the  model  from  zero 
to  maximum  amplitude  rather  than  starting  the  model  with  full  amplitude  at 
the  outset  was  consistent  with  the  slow  closing  procedure  of  the  rotor-to- 
stator  spacing  during  the  actual  test.  In  fact,  attempting  to  start  the 
model  abruptly  almost  always  resulted  in  violently  unsteady  conditions  that 
produced  a model  computation  abort  very  similar  in  nature  to  a surge. 


The  test  procedure  called  for  the  fan  and  P^G  oscillation  to  reach  a 
stable  operating  condition  and  then  the  bypass  discharge  valve  to  be  slowly 
closed  until  surge  was  induced.  Closure  of  the  discharge  valve  produced  an 
increased  blockage  and,  therefore,  increased  the  back  pressure  on  the  fan. 

This  caused  the  fan  operating  point  to  move  up  the  speed  line  characteristic 
either  to  an  operating  condition  that  was  in  equilibrium  with  the  downstream 
conditions  or  to  a point  of  instability.  The  model  performs  precisely  the  same 
function.  The  static  pressure  at  the  discharge  valve  boundary  is  increased 
at  a prescribed  constant  rate  until  aerodynamic  instability  is  encountered. 

The  speed  of  the  model  closure  is  necessarily  faster  than  that  of  the  test 
since  model  operation  is  conducted  on  a microscopic  time  scale  in  comparison 
with  real  time  test  conditions.  It  should  be  noted  that  the  faster  throttling 
rates  used  on  the  model  are  necessitated  by  computation  costs  and  are  justi- 
fied by  the  fact  that  model  rates  of  pressure  increase  of  0.00075  psi  and 
0.00025  psi  per  time  step  have  been  attempted  without  observing  any  large 
dissimilarities  in  results  or  causing  deviation  from  the  steady-state  stage 
characteristics.  A time  step  of  0.0001  sec  was  used  throughout  this  program. 


4 . 6 INSTABILITY  IDENTIFICATION 

It  was  beyond  the  scope  of  this  program  to  determine  if  one  unique 
parameter  would  be  capable  of  identifying  model  instability  and  further, 
would  correlate  with  the  occurrence  of  physical  aerodynamic  surge.  Only 
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individual  flow  parameters  were  examined  as  indications  of  instability  rather 
than  trying  to  form  a unique  combination  of  flow  parameters.  An  effort  to 
predict  instability  using  analytical  expressions  developed  by  other  investi- 
gators (References  16  and  17)  did  not  yield  either  consistent  or  realistic 
results  and  is  omitted  from  discussion. 

Instability  identification  initially  was  completely  based  upon  the 
flow  ratio  and  amplification  function  parameters  developed  during  the  Ref- 
erence 4 study.  The  flow  ratio  parameter,  W2/Wj,  is  a measure  of  the  storage 
of  fluid  within  a volume  since  it  is  the  ratio  of  the  flow  leaving  the  volume 
to  the  flow  entering  the  volume.  A decrease  in  the  flow  ratio  indicates  that 
the  flow  leaving  a blade  row  is  decreasing  at  a faster  rate  than  the  flow 
which  enters  the  blade  row  and  implies  that  the  instability  initiated  in  this 
blade  row.  Upstream  blade  rows  tend  to  show  the  same  response  (V^/W^cl)  and 
the  downstream  blade  rows  show  a flow  ratio  response  greater  than  one  since 
these  blade  rows  are  being  evacuated.  The  amplification  function  is  the  ratio 
of  the  derivative  of  the  volume  averaged  flow  with  respect  to  time  to  the 
derivative  of  the  flow  with  respect  to  time  at  the  system  exit  where  the 
throttling  boundary  condition  is  imposed.  The  amplification  function  is  a 
measure  of  the  rate  at  which  internally  generated  disturbances  grow  with 
respect  to  the  rate  that  the  externally  applied  boundary  condition  is  imposed. 

During  the  course  of  the  study,  three  other  parameters  were  examined  for 
gaining  insight  to  aerodynamic  instability  development. 

1.  The  unsteady  blade  force  which  is  the  component  cf  the  total 
blade  force  associated  with  the  time  rate  of  change  of  angular 
momentum  within  a blade  row. 

2.  The  incidence  angle  which  gives  some  perspective  to  the  stage 
mean-operating  condition  relative  to  the  steady-state  limit  and  the 
degree  of  oscillation  into  the  post-steady-state  surge  region. 

3.  The  mean  overall  operating  point  which  describes  the  mean-overall- 
pressure ratio/mean- inlet-flow  trajectory  on  the  steady-state  fan 
map . 

Discussion  of  the  specific  use  of  these  parameters  and  illustrative 
examples  are  given  in  Section  5.4.2. 
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5.0  RESULTS  OF  DIGITAL  SIMULATIONS  AND  COMPARISONS  WITH  TEST  DATA 


The  results  obtained  from  the  c 1 ean- ini et- flow  throttling  simulations 
and  from  the  planar-wave-distortion  inlet-flow  throttling  simulations  are 
presented  and  discussed  in  detail  in  this  section.  More  specifically,  the 
stage  amplitude  transmission  characteristics,  the  fan  planar  wave  transfer, 
and  stage  initiating  instability  results  have  been  examined  and  are  presented 
here.  Further,  the  simulation  results  obtained  by  including  the  unsteady 
blade  characteristic  effects  are  discussed. 


5 . 1  CLEAN  INLET  FLOW 

Simulation  of  the  clean-inlet  compression-unit  steady-state  performance 
and  establishment  of  the  stability  limit  are  basic  steps  in  the  construction 
of  a dynamic  model.  Simulation  of  the  two-stage  fan  with  clean  inlet  flow 
at  80%  and  100%  corrected  speeds  is  shown  in  Figure  16  relative  to  the 
measured,  non-deteriorated  fan  characteristics.  The  speed-line  character- 
istics pass  through  the  measured  data  points  and  the  stability  limits  are 
coincident  with  the  steady-state  surge  line  as  experimentally  determined. 

The  stability  limit  of  the  steady-state  simulation  at  100%  corrected 
speed  is  defined  by  the  behavior  of  the  blade  row  flow  ratio  W2/W1 . Figure 
17  illustrates  the  blade-row  flow  ratio  (Wy/Wj)  across  both  rotors  of  the 
two-stage  fan  as  a function  of  time.  Applying  the  previously  discussed 
logic  with  respect  to  flow  ratio  behavior  at  surge  (Section  4.6),  it  is 
seen  that  rotor  1 appears  to  be  the  "initiating  stage"  that  produces  the 
fan  instability.  The  rotor  1 exit  flow  is  less  than  the  entrance  flow 
indicating  that  a breakdown  in  the  flow  pumping  capability  of  the  stage  has 
occurred.  Rotor  2 flow  ratio  reacts  to  the  upstream  disturbance  by  reflect- 
ing the  loss  of  flow  through  rotor  1 as  a decrease  in  Wp  relative  to  Wp 
and,  therefore,  a rise  in  W2/W] . 

Figure  18  illustrates  the  behavior  of  the  80%  corrected  speed  steady- 
state  model  instability  characteristic.  The  initiating  stage  is  not  as 
clearly  defined  for  this  surge.  Rotor  1 and  rotor  2 appear  to  go  unstable 
almost  simultaneously.  However,  rotor  2 recovers  while  rotor  1 remains  at 
a low  value  of  Wy/Wp  which  indicates  a more  pronounced  instability. 


5.2  VALIDATION  OF  MODEL  DYNAMIC  CHARACTERISTICS 


5.2.1  Waveform  Fidelity 

A prime  objective  of  this  study  was  to  demonstrate  the  ability  to 
accurately  simulate  the  Planar  Pressure  Pulse  Generator  Program  data  reported 
in  P.eference  3.  To  achieve  that  goal,  certain  test  readings  and  measure- 
ment locations  were  selected  for  comparison  of  model  results  to  test  data 
(See  Figures  2 and  14  and  Tables  3 and  5). 


figure  16 . Clean  Inlet  Flow  Two-Stage  Fan  Map. 


Figure  17.  Airflow  Ratio  Vs.  Time,  10096  N Steady-State  Surg 
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Figure  i8.  Airflow  Ratio  Vs.  Time.  80%  N/\/0,  Steady-State  Surge 


The  test  readings  se  ted  for  simulation  were  originally  chosen  on 
the  basis  of  being  "ope  ng-line"  points.  The  term  "operating-line"  is 
somewhat  nebulous  sin  -tie  test  vehicle  is  a rig  and  not  an  engine.  It 
became  apparent  during  the  modeling  effort  that  a more  convenient  and  cost- 
effective  test  point  could  be  selected.  The  fan  operating  condition  "near 
surge"  was  adopted  as  a more  practical  baseline  point  for  comparison  since 
its  proximity  to  surge  meant  that  considerably  less  computer  time  would  be 
expended  in  throttling  from  the  stationary,  mean  operating-point  up  to  the 
region  of  instability. 

A revision  to  the  waveform  replication  scope  (See  Table  1)  was  made  as 
modeling  at  the  higher  frequencies  progressed.  Analysis  of  the  model 
waveform  character  for  frequencies  greater  than  one-per-rev  (==125  Hz)  indi- 
cated that  the  model  had  limited  capability  to  reproduce  the  higher  harmonic 
components  of  the  fundamental  wave.  At  the  higer  frequencies,  the  waveform 
replication  was  affected  by  the  structure  of  the  wave  pattern  developed  by 
the  Planar  Pressure  Pulse  Generator  and  the  system  boundary  conditions. 

The  major  part  of  the  wave  amplitude  is  associated  with  standing  waves 
rather  than  traveling  waves  in  the  ducting  between  the  P^G  and  the  fan  IGV. 
The  significant  aspect  of  this  type  of  unsteady  flow  is  that  the  wavelength 
at  a given  frequency  is  much  shorter  for  a standing  wave  (Ag)  than  for  a 
traveling  wave  (Ax)  as  shown  by  the  Equation  55. 

Ag  = (1-M)  Ax  (55) 

The  smaller  wavelength  at  a given  frequency  and  the  fixed  "lumped  volume" 
lengths  of  the  model  imply  a reduction  in  spatial  resolution  of  the  waveform, 
especially  for  the  higher  harmonics.  As  a result  of  this  finding  and  in 
recognition  of  the  conclusions  reached  from  the  analysis  of  the  test  data 
indicating  that  stability  losses  were  insignificant  beyond  the  one-per-rev 
frequency,  the  proposed  simulation  matrix  (Table  1)  was  restructed  as 
follows:  1)  All  proposed  simulations  at  frequencies  greater  than  350  Hz 

were  deleted  and,  2)  Additional  emphasis  was  placed  on  low  frequency  analysis 
in  an  effort  to  better  define  the  Planar  Wave  Distortion  Sensitivity  charac- 
teristic. Simulations  were  conducted  at  25  Hz  and  60  Hz,  although  no  test 
data  were  oaken  at  these  frequencies  during  the  experimental  program. 

The  measurement  parameters  that  were  chosen  for  comparison  to  model 
results  have  been  previously  defined  (See  Table  5).  Special  note  should  be 
taken  of  the  fan-face  airflow  parameter  WCA.  This  parameter,  which  is  an 
approximation  to  the  instantaneous  corrected  fan-face  airflow,  is  based  on 
the  instantaneous  physical  airflow  pAV,  the  instantaneous  total  pressure 
Px,  and  the  mean  total  temperature  Tx.  The  mean  temperature  was  used  since 
the  instantaneous  total  temperature  was  not  measured.  This  "approximate" 
corrected  airflow  was  computed  from  the  test  data.  The  identical  set  of 
parameters  were  used  to  calculate  the  same  variable  in  the  model.  During 
the  analysis  of  the  airflow  test  data  during  this  program  it  was  determined 
that  the  original  airflow  data,  as  reported  in  Reference  3,  were  inverted. 
This  inversion  has  been  removed  from  the  experimental  data  used  in  this 
report . 
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The  model  results  that  were  used  for  establishing  the  fidelity  of  the 
model  were  obtained  by  allowing  the  model  to  compute  through  a sufficient 
number  of  cycles  to  insure  that  transients  had  settled  out.  This  procedure 
insured  that  "essentially  stationary"  conditions  were  established  and  pro- 
vided model  results  that  were  representative  of  the  "near  surge"  stationary 
data  point  to  which  the  model  results  would  be  compared.  The  use  of  the 
term  "essentially  stationary"  should  be  noted.  In  some  cases,  both  the 
model  results  and  the  test  data  showed  evidence  of  low-frequency  beat 
frequencies  which  led  to  slight  deviations  from  waveform  stationarity . 

The  results  of  the  waveform  simulations  have  been  summarized  in  terms 
of  amplitude  and  phase  angle  as  a function  of  frequency.  Where  the  model 
amplitudes  (parameter  maximum  minus  parameter  minimum  all  divided  by  the 
average  value  of  the  parameter)  of  the  PTIGV  parameter  did  not  quite  match 
the  test  data,  the  model  values  of  all  parameters  were  scaled  by  the  ratio 
of  the  test  data  PTIGV  value  to  the  model  PTIGV  value.  Examination  of 
Figures  19  and  20  illustrates  some  of  the  aspects  of  model  fidelity.  It  is 
noted  that  at  both  80%  and  100%  corrected  speeds,  there  is  good  agreement 
between  model  and  test  amplitudes  up  to  118  Hz  and  lesser  agreement  above 
this  frequency,  although  the  trends  in  amplitude  are  generally  the  same. 

The  parameter  which  shows  the  greatest  spread  in  values  is  the  P3GPS  param- 
eter. This  probably  is  the  result  of  the  boundary  condition  which  was 
specified  at  the  bellmouth  entrance  and  was  previously  discussed  in  Section 
4.4.3. 

The  phase  angles  of  both  the  model  and  test  parameters  are  measured 
relative  to  the  PTIGV  parameter  and  are  given  in  Figures  21  and  22  for  80% 
and  100%  speeds,  respectively.  Good  agreement  in  the  phase  angle  results 
are  noted  up  to  approximately  80  Hz.  Above  this  frequency,  the  trends  are 
reproduced,  but  there  is  a wide  discrepancy  between  model  and  test  phase 
angles. 

The  phase  angle  plots  of  Figures  21  and  22  may  raise  questions  as  to 
how  waves  propagate,  that  is,  at  lower  frequencies  a parameter  downstream 
of  the  PTIGV  parameter  may  lag  while  at  higher  frequencies  it  may  lead. 
Clarification  of  this  situation  is  difficult  because  both  standing  wave  and 
traveling  wave  components  are  present,  but  it  suffices  to  say  that  the 
phase  angle  relative  to  PTIGV  could  lead  or  lag  the  given  value  by  an 
integral  multiple  of  cycles. 

The  data  comparison  presentations  take  the  following  form:  Eight 

parameters  are  plotted  on  a page  starting  with  the  static  pressure  P3GPS 
and  continue  in  a logical  geometrical  progression  through  the  test  set-up 
concluding  with  the  fan  duct  pressure  at  the  discharge  valve  FDPSD.  Only 
the  unsteady  portion  of  the  waveform  is  shown  since  the  mean  value  of  the 
parameter  has  been  removed.  The  PTIGV  parameter  was  chosen  as  the  reference 
parameter.  The  procedure  for  conducting  waveform  comparisons  begain  by 
overlaying  the  model  waveform  of  PTIGV  on  top  of  the  test  data  PTIGV  wave- 
form, since  all  the  model  waveforms  are  time  correlated  with  each  other  and 
all  the  test  data  waveforms  are  time  correlated  with  each  other.  The  P^G 
model  amplitude  was  chosen  in  a manner  that  caused  as  close  a match  as 
practical  between  the  PTIGV  model  amplitude  and  the  PTIGV  test  data  ampli- 
tude. Hence,  the  degree  of  waveform  similitude  between  the  model  results 
and  the  test  data  can  be  assessed  for  the  other  parameters. 
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5.2. 1.1 


Simulation  of  Sdg.  71,  1007'.  N//o,  42  Hz 


Figure  23  illustrates  the  comparison  between  the  model  resuits  and  the 
test  data  for  100%  corrected  speed  and  a Planar  Pressure  Pulse  Generator 
frequency  of  42  Hz.  The  model  simulates  the  amplitude  and  phase  of  all 
variables  throughout  the  system.  The  model  static  pressure  forward  of  the 
P^G  (P3GPS)  shows  a higher  amplitude  than  test  results.  This  variation  is 
attributed  to  the  forward  boundary  condition  (P-p  and  Tp  held  constant)  as 
previously  discussed  (See  Section  4.4.3),  but  is  not  considered  to  have 
significant  impact  on  the  overall  accuracy  of  the  simulation.  Both  the 
model  and  the  test  data  PS14  parameter  exhibit  second  harmonic  content 
although  the  test  data  exhibits  it  to  a lesser  extent.  It  should  be  noted 
that  the  model  is  being  slowly  throttled  and  that  the  amplitude  of  the  PS14 
and  FDPSC  are  both  being  attenuated  by  the  change  in  the  operating  point 
boundary  condition.  Since  it  was  not  always  possible  to  precisely  reproduce 
the  test  mean  operating  point,  some  deviation  can  be  anticipated  as  a 
result  of  this  operating  point  d ssimilarity.  The  model  "approximate" 
airflow  matches  the  test  data  "approximate"  airflow. 


5. 2. 1.2  Simulation  of  Rdg.  82,  80%  N//9,  42  Hz 

Figure  24  illustrates  the  80%  corrected  speed,  42  Hz  validation  case. 
Results  are  similar  to  the  100%  speed,  42  Hz  case  with  the  exception  of  the 
presence  of  a strong  second  harmonic  contribution  to  both  the  model  and 
test  "approximate"  corrected  airflow  and  pressure  waveforms  downstream  of  the 
fan.  The  amplitudes  of  all  model  parameters  show  agreement  with  the  test 
data.  The  amplitude  and  phase  of  the  model  second  harmonic  contribution  is 
not  in  phase  with  the  test  values,  but  the  shape  of  these  waveforms  are  not 
entirely  stationary.  A model  simulation  of  longer  duration  would  probably 
produce  l-etter  agreement  between  waveforms. 


5. 2. 1.3  S i mujation  o_f_  _Rdg_.  69,  _100%  N//e,  80_Hz 

Figure  25  shows  the  validation  curves  for  this  case.  The  standing 
wave  amplitudes  and  phase  are  illustrated  by  the  static  pressures  in  the 
duct  between  the  P^G  and  the  fan  (AXPS2  and  AXPS3) . The  AXPS2  measurement 
appears  to  be  near  a node  in  the  pressure  wave  since  the  amplitudes  are 
small  and  the  waveform  is  erratic.  AXPS3  is  a pressure  measurement  approxi- 
mately 15  inches  downstream  of  AXPS2  and  shows  large  pressure  oscillations 
of  uniform  amplitude  which  are  indicative  of  being  closer  to  an  antinode. 

The  model  results  replicate  the  phase  and  amplitude  of  these  waveforms. 

The  test  data  "approximate"  corrected  airflow  was  not  available  for  this 
reading.  There  is  only  a slight  deviation  between  model  and  test  amplitudes 
of  the  downstream  pressures  while  the  phase  relationships  remain  well 
correlated . 
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Figure  25.  Simulation  of  lidy.  69,  100%  N/v^  80  Hz 
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5 . 2 . 1 . 4  Simulation  of  Rdg.  77,  80%  N//e,  75  Hz 

Figure  26  illustrates  the  test  and  model  results  for  this  case.  The 
duplication  of  phase  and  amplitude  of  pressures  throughout  the  simulation 
is  shown  to  be  excellent.  The  model  corrected  airflow  phase  and  magnitude 
is  shown  to  deviate  somewhat  from  the  test  data.  No  apparent  explanation 
exists  for  this  departure,  but  the  experimental  nature  of  the  airflow 
measurement  device  must  be  considered  when  assessing  the  quality  of  this 
test-data  measurement. 


5 . 2 . 1 . 5 Simulation  of  Rdg.  64,  100%  N/Sd,  118  Hz 

Figure  27  represents  the  simulation  of  this  case.  The  118  Hz  frequency 
test  data  contained  large  harmonic  components  that  distort  the  fundamental 
waveform.  The  model  does  not  simulate  these  higher  frequency  components. 
Model  and  test  amplitudes  compare  reasonably  well  over  the  entire  system, 
but  phase  shifts  are  apparent  throughout.  As  frequencies  increase  and 
wavelengths  decrease,  the  positions  of  the  instrumentation  relative  to  each 
other  (model  versus  test)  and  to  the  boundary  conditions  become  important 
to  the  precision  of  the  phase  comparisons,  especially  when  dealing  with 
standing  wave  components. 

This  case  ^lso  provides  the  base  point  for  two  other  simulations. 

Model  operation  at  two  different  pressure  fluctuation  amplitudes  was  carried 
out  for  the  purposes  of  establishing  distortion  sensitivities  as  a function 
of  amplitude.  The  results  of  the  sensitivity  variation  with  amplitude  will 
be  discussed  later. 

5.2. 1.6  Simulation  of  Rdg.  61,  80%  N//6,  118  Hz 

Figure  28  illustrates  the  model  versus  test  comparisons  for  this 
reading.  The  comments  addressed  to  the  previous  118  Hz  point  can  also  be 
similarly  applied  to  the  80%  corrected  speed  case.  However,  the  contribution 
of  the  second  harmonic  to  the  overall  test  waveform  is  somewhat  less  than 
in  the  previous  case  and  the  phase  relationships  are  improved.  An  exception 
to  the  overall  agreement  is  the  P3GPS  model  results  with  its  larger  ampli- 
tude and  little  evidence  of  higher  harmonic  frequencies  as  compared  to  the 
test  data.  Other  variables  (AXPS2  and  PTIGV)  show  similar  effects,  but 
with  a reduced  harmonic  amplitude. 
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5. 2. 1.7  Simulation  of  Rdg.  41,  100%  N//8,  220  Hz  and  Rdg.  57, 

80%  ¥/7o",  240  Hz 

Simulations  of  Rdgs.  41  and  57  are  shown  in  Figures  29  and  30,  respec- 
tively. The  model  results  show  responses  similar  to  the  test  data  responses 
for  both  cases  at  these  frequencies.  The  fundamental  wave  shapes  and 
amplitudes  are  preserved,  but  the  harmonics  of  the  test  data  are  not  seen  in 
the  model  results.  Phase  relationships  are  not  maintained  as  well  as  with 
the  lower  frequency  cases. 

The  amplitude  and  phase  of  the  predicted  airflow  for  the  Rdg.  41  case 
are  in  agreement  with  the  amplitude  and  phase  of  the  measured  airflow, 
despite  the  non-uniformity  of  the  measured  airflow  waveform.  Phase  agree- 
ment is  to  be  expected,  since  the  measurement  planes  for  both  the  model  and 
the  test  data  airflows  are  coincident. 

The  Rdg.  57,  model  total-pressure  amplitude  at  the  IGV  station  does 
not  match  the  measured  PTIGV  value  because  it  was  necessary  to  hold  the 
model  amplitude  at  a lower  value  to  insure  stable  model  operation  at  this 
condition.  This  apparent  sensitivity  of  the  model  is  attributed  to  the 
extremely  large  variation  in  "approximate"  corrected  airflow  at  the  fan 
inlet  (WCA  amplitude  =30%).  Measured  airflows  were  not  available  at  this 
condition  for  comparison. 


5. 2. 1.8  Simulation  of  Rdg.  43,  100%  N//6,  350  Hz  and  Rdg.  53, 
80%  N/Ve  ,'~350  Hz 


Simulations  of  Rdgs.  43  and  53  are  shown  in  Figures  31  and  32,  respec- 
tively. The  most  significant  feature  of  these  comparisons  is  the  apparent 
attenuation  of  model  static  pressures  as  the  distance  from  the  Planar  Pres- 
sure Pulse  Generator  increases.  At  this  frequency  for  both  the  100%  and 
80%  conducted  speeds,  the  model  P3GPS  parameter  shows  an  attenuation  of 
pulse  amplitude  relative  to  test  data,  while  at  lower  frequencies  the  model 
P3GPS  parameter  had  tended  to  exhibit  amplitudes  greater  than  these  of  the 
test  data.  Attenuation  of  the  amplitudes  is  thought  to  be  a function  of 
the  model  transmission  characteristics  at  this  fundamental  frequency,  in 
the  same  manner  as  the  higher  harmonics  of  the  lower  frequencies  appeared 
to  be  attenuated. 


5.2.2  Sample  Plotted  Results 

The  overall  analysis  of  each  model  case  was  facilitated  by  the  use  of 
machine  plots  of  pertinent  aerodynamic  variables  displayed  as  a function  of 
time.  A set  of  these  plots  are  shown  in  Figures  33  through  39  for  the  Rdg. 
82,  100%  corrected  speed,  42  Hz  case  simulation. 

Figure  33  shows  the  variation  of  physical  airflow  entering  each  rotor 
as  a function  of  time.  The  phase  differences  in  flow  are  suggestive  of 
stage  mismatch  in  flow  handling.  Figure  34  illustrates  the  variation  of 
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the  incidence  angle  as  a function  time  for  both  rotor  inlets.  As  indicated 
in  Section  A. 4. 2,  the  incidence  angle  is  functionally  related  to  the  relative 
total-pressure  loss  coefficient  and  the  deviation  angle  parameters.  Figures 
35  and  3b  respectively,  illustrate  the  nonlinearities  introduced  by  these 
parameters.  Specifically,  the  rotor  1 incidence  angle  (Figure  34)  is  exper- 
iencing low  level  oscillations  which  occur  on  the  flat  part  of  the  loss 
coefficient-incidence  angle  curve  and  thus,  give  rise  to  very  low  levels  of 
loss  coefficient  oscillations.  On  the  other  hand,  the  rotor  2 incidence 
angle  is  going  through  a significant  range  which  happens  to  fall  on  a rising 
portion  of  the  loss  coefficient-incidence  angle  curve  and  thus  introduces 
large  variations  in  loss  coefficient.  Figures  37,  38,  and  39  illustrate  the 
time-dependent  variation  in  rotor  total-pressure  ratio,  static-pressure 
ratio,  and  total-temperature  ratio,  respectively.  The  sets  of  plots  for  all 
model  simulation  baseline  points  are  presented  in  Appendix  C. 

5 . 3 FAN  I’LANAR  WAVE  DISTORT  ION  TRANSFER 

The  detailed  fan  model  results  permit  the  investigation  of  unsteady 
total-pressure  wave  propagation  through  the  fan  test  vehicle.  The  amplitude 
transmission  ratio  parameter  established  for  analyzing  these  results  is 
defined  as  the  nondimens ional  total-pressure  variation  ( (max-min) /avg)  at 
the  downstream  station  divided  by  the  nondimensional  total  pressure  varia- 
tion ( (max-min) /avg)  at  the  upstream  station.  The  stations  chosen  for  this 
analysis  were:  The  1GV  entrance,  rotor  1 entrance,  rotor  2 entrance  and 

plane  14.  The  fan  sections  which  were  analyzed  and  the  associated  ampli- 
tude transmission-ratio  designations  are  tabulated  below. 

IGV  - PTri/Ptigv 

Stage  1 - PTR2/pTri 

Stage  2 - pTi4/pTR2 

Overall  Fan  - pTi4/pTlGV 

The  amplitude  transmission  ratios  for  baseline  and  near  instability  points 
were  determined  for  each  model  test  case. 

Figure  40  illustrates  the  computed  amplitude  transmission  character- 
istics for  the  100%  corrected  speed,  42  Hz  frequency  simulations  of  Rdg. 

71.  As  shown,  the  stage  2 characteristic  shows  a marked  decrease  when 
throttling  from  the  baseline  points  towards  instability  and  is  the  major 
contributor  to  the  attenuation  exhibited  by  the  overall  transfer  character- 
istic pTl4/pTiGV"  The  detailed  amplitude  transmission  characteristics  for 
each  model  test  case  are  given  in  Appendix  D.  riots  such  as  these  graphi- 
cally provide  information  on  which  parts  of  a compression  component  are 
amplifying  or  attenuating  unsteady  disturbances. 
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Figures  41  through  43  illustrate  the  internal  transfer  characteristics 
as  .1  function  of  frequency  for  1 00%  corrected  speed  operation  for  the  iC.V, 
stage  1,  and  stage  2,  respectively.  Values  of  amplitude  ratio  greater  than 
one  occur  across  the  IGV.  The  existence  of  this  amplification  is  not 
readily  understood;  however,  other  model  and  test  data  for  other  compres- 
sion components  tend  to  support  this  result.  Stage  1 (Figure  42)  is  shown 
to  generally  attenuate  the  unsteady  waves  at  frequencies  below  100  Hz  and 
at  ISO  Hz.  The  118  Hz  and  220  Hz  points  indicate  that  amplification  of 
entering  flow  disturbances  would  occur  in  this  frequency  range.  The  stage 
2 amplitude  ratio  (Figure  43)  varies  significantly  across  the  frequency 
range.  Transmission  ratios  are  1.0  or  less  in  all  cases  with  the  25  Hz 
case  falling  as  low  as  0.18. 

The  80%  corrected  speed  IGV,  stage  1,  and  stage  2 amplitude  ratios  as 
a function  or  frequency  are  shown  in  Figures  44  through  4G.  The  IGV  ampli- 
tude ratio  is  again  equal  to  1.0  or  greater  in  most  cases  with  the  near 
instability  data  showing  an  increase  in  transfer  in  every  case  relative  to  the 
baseline  point.  Stage  1 attenuation  increases  with  increasing  frequency  and 
the  near  surge  data  demonstrates  somewhat  lower  transfer  coefficients  than 
the  operating  line  data.  The  stage  2 transfer  shows  that  attenuation  exists 
across  this  stage  for  ali  frequencies  except  at  118  hz. 

Near  instability  conditions  as  compared  to  baseLine  conditions,  arc* 
not  widely  different  operational  points.  The  80%  corrected  speed  transfer 
results  are  the  more  definitive  set  since  these  points  were  generally 
throttled  further  in  terms  of  flow  from  the  baseline  point  to  near  insta- 
bility. 

The  overall  planar  wave  distortion  transmission  properties  of  the  fan 
are  defined  by  the  ratio  of  the  nondimensionalized  amplitude  [ (max-min) /avgl 
of  the  exiting  (plane  14)  total-pressure  wave  to  the  nondimensionalized  ampli- 
tude of  the  entering  (IGV  plane)  total-pressure  wave.  Fan-planar-wave 
transfer  was  determined  from  test  data  for  both  the  hub-  and  tip  flow  of 
the  fan.  The  hub  transfer  was  based  on  data  taken  at  radial  immersions 
corresponding  to  10%  and  30%  of  the  flow  area  and  the  tip  transfer  was 
based  on  data  taken  at  radial  immersions  of  70%  and  90%  of  the  flow  area. 

The  test  data  indicated  a significant  change  in  transfer  properties  witli 
respect  to  radial  location  and  since  the  model  calculations  are  based  on 
pitch  line  characteristics,  comparison  of  model  and  test  results  must 
include  both  hub  and  tip  transfer  coefficients  from  the  test  data. 

The  transfer  characteristics  as  a function  of  frequency  for  100% 
corrected  speed  are  illustrated  in  Figure  47.  The  genera]  behavior  of  the 
model  results  matches  the  test  data  and  falls  between  the  hub  and  tip 
results  as  might  be  anticipated.  The  model  results  indicate  large  variations 
in  transfer  properties  at  low  frequencies.  The  transfer  coefficients  at  25 
Hz  and  60  Hz  are  very  small  but  at  42  Hz  an  increased  transfer  of  approxi- 
mately 300%  above  both  the  25  Hz  and  60  Hz  cases  is  noted. 
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Figure  48  illustrates  the  80%  corrected  speed  overall  transfer  coef- 
ficient comparisons.  The  transfer  "spike"  at  118  Hz  shown  by  the  test  data 
is  the  predominate  feature  of  the  characteristic  and  is  reproduced  accurately 
by  the  model.  Although  these  results  suggest  that  the  fan  may  act  as  an 
amplifier  at  118  Hz,  this  interpretation  of  the  results  would  probably  be 
misleading.  Rather,  it  is  suspected  that  the  results  arise  from  the  presence 
of  standing  waves  in  the  ducting  upstream  and/or  downstream  of  the  fan. 

The  planar  wave  distortion  transfer  results  are  summarized  in  Table  6. 

It  is  noticed  that  the  model  overall  transfer  values  lie  between  the  hub 
and  tip  test  values  for  all  cases  but  one.  These  results  were  expected 
since  the  model  is  a pitch  line  model. 

Table  6.  Summary  of  Model  and  Test  Planar  Wave  Distortion  Transfer  Results. 


Planar  Wave  Distortion  Transfer 


Simulated 

Model 

Test 

Reading  No. 

Frequency 

IGV 

Stage  1 

Stage  2 

Overall 

Hub 

He 

100%  Speed 

- 

25 

1.000 

0.702 

0.177 

0.125 

- 

- 

71 

42 

1.074 

0.  792 

0.706 

0.601 

0.19 

0.45 

- 

60 

1.006 

0.  560 

0.379 

0.180 

- 

- 

69 

80 

1.148 

0.767 

0.766 

0.676 

0.36 

1.12 

64 

118 

1.119 

0.842 

0.907 

0.907 

0.13 

1.03 

41 

220 

0.923 

1.433 

0.373 

0.494 

0.18 

0.91 

43 

350 

0.970 

0.  838 

0.622 

0.505 

0.17 

0.72 

80%  Speed 

82 

42 

1.002 

0.954 

0.772 

0.739 

0.22 

0.80 

77 

75 

0.  967 

0.  890 

0.610 

0.520 

0.15 

0.  92 

61 

118 

1. 148 

0.765 

1.453 

1.280 

1.21 

2.14 

57 

240 

1.  423 

0.619 

0.888 

0.783 

0.47 

1.02 

66 

350 

1.102 

0.712 

1.048 

0.823 

0. 17 

1.08 

The  fact. 

that  the 

model  predicts  the 

abnormally  large  transfer 

at  80% 

corrected  speed,  118  Hz 

and  that 

the  model 

results 

follow  the 

trends 

of  the 

test  data  with 

the  model 

values 

falling  between  hub 

and  tip  test  results, 

is  clear  evidence  of  the  validity  and  capability  of  this  model  to  properly 
handle  the  interaction  of  planar  waves  with  a compression  component. 
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5.4  INSTABILITY  ANALYSIS 


The  results  obtained  by  throttling  the  model  to  instability  are  presented 
In  this  section.  Depending  upon  corrected  speed  and  planar  wave  frequency, 
it  has  been  possible  to  characterize  four  distinctly  different  types  of 
model  instability  which  appear  to  have  counterparts  during  actual  vehicle 
operation.  These  types  of  instability,  planar  wave  distortion  sensitivity, 
and  stage  initiating  instability  results  are  presented  in  the  following 
paragraphs . 


5.4.1  General  Observat ions 

The  similarity  between  model  behavior  when  aerodynamic  instabilities 
occur  and  actual  fan  surge  is  the  most  significant  observation  of  the 
entire  instability  point  analysis.  Model  response  parallels  the  trends 
defined  by  test  data  in  several  areas.  The  "hard"  surges  associated  with 
high  speed  fan  operation  and  the  "softer"  surges  seen  during  medium  and  low 
speed  operation  are  reproduced  by  instabilities  which  develop  within  the 
model.  Non-stationary  events  which  may  be  analogous  to  rotating  stalls  are 
observed  in  the  model  calculations.  Finally,  migration  of  the  fan  model 
mean  operating  point,  especially  near  instability  at.  low  speeds,  shows 
close  qualitative  agreement  with  test  data. 


5.4.2  Discussion  of  Instability  Point  Data 

The  analysis  of  the  model  instability  point  is  based  upon  the  param- 
eters discussed  in  Section  4.6  Instability  Identification.  The  instability 
encountered  when  throttling  the  clean-inlet  model  was  discussed  under  Clean 
Inlet  Flow  results  in  Section  5.1.  The  blade-row  flow  ratio  technique  for 
determining  the  stage  initiating  instability  and  the  point  of  instability 
was  illustrated  for  the  80%  and  100%  speed  corrected-bpeed  lines.  This 
instability  has  been  termed  a Type  1 Instability.  Analysis  and  illus- 
tration of  each  of  the  remaining  types  of  instability  are  presented  below. 


Type  2 Instability 

This  type  of  instability  occurred  at  100%  corrected  speed  for  planar 
wave  frequencies  in  the  range  of  42  Hz  to  220  Hz.  All  of  these  instabili- 
ties developed  very  abruptly  in  the  model.  The  migration  of  the  mean 
operating  point  on  the  steady-state  fan  map  was  nominally  well  behaved 
until  the  point  of  model  instability  was  reached.  Model  computational 
instability  occurred  within  a few  time  steps  after  aerodynamic  instability 
occurred  for  the  majority  of  these  cases. 

The  Type  2 Instability  is  illustrated  by  the  80  Hz  simulation  at  100% 
corrected  speed.  The  stage  initiating  instability  was  derived  from  the 
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Wi/Wi  relationship  and  was  substantiated  by  the  response  of  the  "unsteady” 
blade  force.  Figure  49  and  Figure  50  illustrate  the  behavior  of  V^/W^  and 
the  blade  force  as  a function  of  time,  respectively.  The  flow  ratios  across 
both  rotors  are  periodically  higher  than  the  values  at  steady-state  surge 
without  producing  instability.  The  periodicity  of  the  imposed  unsteady 
fluctuations  permits  these  high  flow  ratios  to  exist  without  surge  until  a 
large  deviation  occurs  in  rotor  2 as  seen  in  Figure  49.  The  unsteady  blade 
force  term  (Figure  50)  shows  similar  behavior  in  rotor  2 where  an  apparently 
instantaneous  breakdown  is  observed  as  a sharp  rise  and  marks  rotor  2 as  the 
initiating  stage. 


Type  3 Instability 

This  type  of  instability  occurred  at  80%  corrected  speed  for  planar 
wave  frequencies  in  the  range  of  42  Hz  to  240  Hz.  This  class  of  instability 
differs  considerably  from  the  100%  corrected  speed  instability  character- 
istics in  the  same  frequency  range.  None  of  these  cases  experienced  a 
numerical  instability  that  would  cause  the  solution  to  abort.  Large  fluc- 
tuations occur  in  the  internal  aerodynamics  of  the  fan,  but  the  mean  per- 
formance is  well  behaved  until  the  fan  is  throttled  higher  on  its  speed 
line.  The  point  of  instability  for  this  class  of  instability  is  defined  by 
the  last  position  of  the  mean  operating  conditions  prior  to  initiation  of 
erratic  mean  behavior. 

A Type  3 Instability  is  illustrated  by  Figure  51  for  the  240  Hz,  80% 
speed  case.  The  operating  condition  at  the  abrupt  change  in  slope  of  the 
speed  line  is  defined  as  the  instability  point  for  the  Type  3 Instability. 

The  initiating  stage  for  this  type  of  instability  may  be  best  defined 
by  the  amplification  function.  This  parameter  is  defined  as  the  ratio  of 
the  time  derivative  of  the  volume  averaged  flow  in  a given  rotor  to  the 
ratio  of  the  time  derivative  of  flow  in  the  exit  volume  where  the  throttling 
process  is  controlled  (AMPF  = DW(x)rotor/DW(x)exit) • Figure  52  represents 
the  behavior  of  this  function  for  the  period  of  time  associated  with  the 
development  of  the  Type  3 Instability.  It  is  clearly  shown  that  the  rotor 
2 amplification  function  exhibits  large  excursions  relative  to  the  rotor  1 
amplification  function.  On  the  basis  of  these  observations,  rotor  2 is 
determined  to  be  the  stage  initiating  instability  for  this  case. 


Type  4_ Instability 

This  type  of  instability  occurs  at  350  Hz  at  both  80%  and  100%  cor- 
rected speeds.  The  point  of  model  instability  for  both  of  these  cases  was 
determined  to  be  very  close  to  the  steady-state  clean-inlet-surge  line. 

The  behavior  of  these  points  was  such  that  a clear  determination  of  the 
initiation  of  surge  could  not  be  made.  Occasional  disturbances  which  may 
be  physically  analagous  to  rotating  stalls  were  observed  in  the  data  as 
shown  in  Figure  53  which  represents  the  near  instability  portion  of  the 
100%  corrected  speed,  350  Hz  simulation. 
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Figure  lo.  Unsteady  Blade  Force  Vs.  Time,  IOO%  N/\^  80  Hr, , 


The  path  of  the  mean  operating  point  for  this  case  is  shown  in  Figure 
54.  The  speed  line  described  by  the  mean  operating  point  "rolls-over"  and 
moves  down  the  positive  slope  portion  of  the  characteristic  for  some  distance 
prior  to  model  -calculat ion  abort.  The  point  chosen  to  represent  instability 
was  the  zero  slope  position.  This  choice  is  consistent  with  the  interpreta- 
tion of  steady-flow  fan  stability,  speed  line  slope,  and  with  previous  test 
data  analyses.  Figure  54  also  shows  the  path  described  by  the  fan  mean 
operating  point  based  upon  analysis  of  the  high-response  test  data.  The 
mean  operating  point  rotor  described  by  the  average  of  the  dynamic  test 
data  may  not  be  quantitatively  accurate,  but  its  trend  appears  to  be  valid. 
Hence,  the  trend  of  model  results  agrees  with  the  trend  of  the  dynamic-data 
calculation. 

The  initiating  stage  for  this  class  of  instability  was  not  determined. 
Figure  55  illustrates  the  growth  of  the  unsteady  force  term  as  the  mean 
point  moves  into  the  region  of  instability  but,  it  does  not  yield  a defini- 
tive evaluation  of  the  initiating  stage. 

The  results  of  characterizing  the  types  of  instability  are  summarized 
in  Figures  56  and  57  where  the  pressure  ratio-corrected  flow  relationships 
and  flow  parameters  are  presented,  respectively  for  each  type  of  instability. 

The  Type  1 Instability  is  characterized  by  a relatively  slow  bending 
over  of  the  speed  line  for  clean-inlet-flow  conditions.  The  flow  ratio 
W^/Wj  is  the  parameter  most  useful  for  determining  the  stage  which  initiates 

instability. 

The  Type  2 Instability  is  characterized  by  the  abrupt  development  of 
instability  and  is  associated  with  planar  wave  frequencies  in  the  range  of 
42  Hz  to  220  Hz  at  100%  corrected  speed.  Although  the  flow  ratio  is 

very  useful  for  determining  the  stage  which  initiates  instability,  the 
umpl if ication  function  also  provides  similar  information.  The  data  shown  on 
Figure  57  illustrate  the  abrupt  development  of  the  instability. 

The  Type  3 Instability  is  characterized  by  erratic  behavior  of  the 
mean  operating  point  in  terms  of  the  pressure  ratio  - corrected  flow  relation 
ship  and  is  associated  with  planar  wave  frequencies  in  the  range  of  42  Hz 
to  240  Hz  at  80%  corrected  speed.  The  amplification  function  is  the  most 
useful  paiameter  for  determining  the  stage  where  an  instability  originates 
by  the  large  magnitude  of  the  erratic  fluctuations  the  initiating  stage 
experiences . 

The  Type  4 Instability  is  characterized  by  the  rapid  bending  over  of 
the  speed  line  at  350  Hz  for  both  80%  and  100%  corrected  speeds.  None  of 
the  parameters  that  were  being  examined  during  the  course  of  this  study 
were  useful  in  determining  the  stage  where  an  instability  of  this  type 
was  initiated.  However,  behavior  was  noted  in  the  amplification  parameter 
which  might  be  associated  with  rotating  stall,  that  is,  anti-phasing  between 
stages  was  noted  during  the  first  few  cycles  as  shown  on  Figure  57. 
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5.4.3  Surge  Pressure  Ratio  Sensitivity 

The  loss  in  surge  pressure  ratio  per  unit  amplitude  of  inflow  planar 
wave  distortion  is  defined  as  surge  pressure  ratio  sensitivity.  Loss  in 
surge  pressure  ratio  (APRS)  is  defined  as  the  incremental  change  in  the 
pressure  ratio  at  surge  as  measured  from  the  clean-inlet  surge  line  at 
constant  corrected  airflow  to  the  surge  line  with  distorted  inflow  divided  by 
the  clean-inlet-surge-line  pressure  ratio.  The  planar  wave  distortion  ampli- 
tude is  defined  by  the  following  relationship: 
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Therefore,  the  distortion  sensitivity  to  planar  waves  is  defined  to  be 

DS  = APRS/[APT/(2PT)IGV]  (57) 

The  definition  of  loss  of  surge  pressure  ratio  with  planar  wave  distortion 
requires  additional  constraints  in  order  to  specify  the  surge  pressure 
ratio  with  inflow  planar  waves  since  the  instantaneous  fan  operating  point 
is  not  fixed.  The  loss  in  surge  pressure  ratio  was  determined  using  the 
mean  operating  point  of  the  fan  just  prior  to  instability  as  determined 
from  steady-state  instrumentation.  The  analogous  procedure  in  the  model  is 
to  time  average  th  : instantaneous  pressure  ratio  and  flow  values  for  each 
cycle.  As  the  model  is  throttled,  a locus  of  the  average  values  is  scribed 
on  the  steady-state  fan  map  and  represents  the  model  version  of  the  progres- 
sion of  the  fan  mean  operating  point.  The  mean  point  at  instability  was 
selected  based  upon  the  instability  analyses  that  were  conducted  according  to 
the  discussion  of  Section  5.4.2.  Comparison  of  the  model  100%  corrected- 
speed  planar-wave  distortion  sensitivities  with  the  test  results  is  shown  in 
Figure  58.  The  trend  of  model  calculated  sensitivity  is  shown  to  be  similar 
to  the  test  results.  Two  additional  model  points  are  shown  for  which  test 
data  does  not  exist  (25  Hz  and  60  Hz).  These  cases  were  added  to  the  simu- 
lation matrix  in  order  to  develop  the  sensitivity  characteristic  at  lower 
frequencies.  The  shape  of  the  model  sensitivity  characteristic  in  this 
region  suggests  that  this  function  is  more  non-linear  with  frequency  than  had 
been  concluded  from  the  earlier  test  results. 

A consideration  not  previously  discussed  with  respect  to  sensitivity 
is  the  effect  of  deterioration  or  the  test  sensitivity  results.  The  poten- 
tial variation  in  the  42  Hz  sensitivity  is  also  shown  on  Figure  58  if  the 
post-test  surge  line  is  used  instead  of  the  pre-test  surge  line.  Since 
this  surge  came  near  the  end  of  the  test,  the  inclusion  of  the  deterioration 
effect  should  be  more  correct. 
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The  effect  of  varying  amplitude  on  dynamic  distortion  sensitivity  at  a 
given  frequency  is  shown  also  in  Figure  58  at  118  Hz.  Total-pressure 
amplitudes  of  8.20,  9.66,  and  11.40  percent  were  evaluated  for  point  of  in- 
stability and  dynamic  distortion  sensitivity.  Sensitivities  of  0.634,  0.673, 
and  0.640,  respectively,  were  computed  for  these  total-pressure  amplitudes. 
When  normal  data  scatter  and  the  effect  of  the  small  magnitude  of  the  surge 
pressure  ratio  loss  in  the  computation  are  taken  into  account,  the  variability 
of  test  data  sensitivity  does  not  nave  a strong  relationship  with  planar  wave 
amplitude. 

The  overall  levels  of  model  and  test  sensitivity  are  in  good  agreement 
with  the  exception  of  the  frequency  range  between  100  and  300  Hz.  It  is 
possible  that  part  of  the  variation  seen  in  this  frequency  range  is  due  to  the 
differences  between  the  extrapolated  model-blade-row  characteristics  in  the 
post-surge  region  and  the  actual  fan  characteristics.  In  addition,  the 
unsteady  flow  effects  on  blade  characteristics  become  increasingly  important 
in  this  region  and  the  model  analyses  have  not  taken  this  factor  into  account. 
The  decrease  in  measured  sensitivity  with  increasing  frequency  until  the 
effect  of  the  unsteady  flow  on  stability  is  negligible  is  reproduced  by  the 
model  results. 

The  80%  speed  planar  wave  distortion  sensitivity  model  and  test  results 
are  shown  in  Figure  59.  Model  sensitivities  simulate  test  results  across  the 
frequency  range  with  the  possible  exception  of  the  42  Hz  case.  As  with  the 
1 007  speed,  42  Hz  case,  the  deterioration  effect  may  be  taken  into  account  as 
shown  in  Figure  59.  The  .80%  model  sensitivity  characteristic  shape  has  not 
been  explored  at  other  additional  low  frequencies  as  was  done  for  100% 
corrected  speed,  but  it  is  assumed  that  its  behavior  would  be  similar  to  the 
100%  corrected  speed  sensitivity  characteristic  at  frequencies  below  42  Hz. 

The  amplitudes  used  to  determine  the  planar  wave  sensitivities  and  the 
sensitivities  are  tabulated  in  Table  7 for  both  the  model  results  and  the  test 
data. 


5.4.4  Sample  Plotted  Results 

The  instability  analyses  were  based  largely  upon  a selected  group  of 
parameters  which  exhibited  a strong  functional  relationship  to  the  breakdown  of 
stable  fan  operation.  These  parameters  were  the  unsteady  blade  force,  the 
amplification  function,  the  tangent  of  the  incidence  angle,  and  the  flow  ratio. 
The  unsteady  blade  force  behavior  is  shown  in  Figure  60  for  118  Hz,  100% 
corrected  speed  as  the  model  enters  the  region  of  aerodynamic  instability. 

The  unsteady  force  term  is  a function  of  the  blade  row  volume  dynamics  and  is 
sensitive  to  destabilizing  forces.  An  example  of  the  behavior  of  the  amplifica- 
tion function  is  shown  in  Figure  61  for  the  100%  speed,  118  Hz  case.  Figure 
62  illustrates  the  instantaneous  behavior  of  the  tangent  of  the  incidence 
angle  at  the  entrance  to  each  rotor.  Observation  of  this  parameter  relative 
to  the  location  of  the  steady-state  surge  incidence  angle  gave  some  insight 
into  the  blade  row  characteristic  fluctuations  relative  to  steady-state  opera- 
tion. Each  of  these  parameters  are  shown  for  each  surge  in  Appendix  E as  an 
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Table  7.  Model  and  Test  Planar  Wave  Distortion  Sensitivities. 
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Simulated 
Reading  No. 

Frequency 

Planar  Wave 
Amplitude 

Planar  Wave  Distortion 
Sensitivity 

Model 

Test 

Model 

Test 

1002  N/_/e 

- 

25 

7.20 

- 

0.803 

- 

71 

42 

7.40 

6.88 

0.432 

0.74/0.651* 

- 

bO 

11.82 

- 

0.333 

- 

69 

eo 

10.80 

10.67 

0.682 

0.59 

64 

118 

8.  20 

8.79 

0.634 

0.14 

64 

118 

9.66 

- 

0.673 

- 

65/66 

118 

11.40 

15.39 

0.640 

0.39 

41 

220 

7.07 

6.06 

0.565 

-0.28 

43 

350 

8.70 

10.48 

0.1h2 

-0.11 

80%  N//e 

82 

42 

6.70 

5.84 

0. 164 

0.480/0. 362* 

77 

7 5 

6.05 

7.09 

0.152 

-0.02 

61 

118 

5.99 

5.11 

0.230 

0.14 

57 

240 

6.36 

9.66 

0.044 

0.07 

53 

350 

6.05 

6.41 

-0.139 

-0.35 

*Potential 

deter Iorat 

ion  taken 

into  acount 
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Figure  6l.  Amplification  Function  Vs 
N/v/S’,  118  Hz. 


TAN(i) 


aid  to  the  interpretation  of  the  model  instability.  The  flow  ratio  parameter 
(W2/W^)  was  not  always  available  in  plotted  form,  but  was  available  in  the 
computer  listings  for  use  in  the  analyses. 


5.4.5  Stage  Initiating  Instability 

The  stage  initiating  aerodynamic  instability  was  determined  by  the  criteria 
set  forth  in  Section  5.4.2  on  instability  identification  and  with  the  aid  of 
the  plotted  parameters  described  in  Section  5.4.4.  Based  upon  analysis  of  each 
of  these  parameters,  the  stage  initiating  instability  was  determined  where 
positive  identification  could  be  made  and  the  results  are  presented  in  Table  8. 
Behavior  of  the  individual  stages  fall  into  a general  pattern.  Stage  two 
activity  level  was  usually  high  in  comparison  to  stage  one  and  as  the  fan  was 
throttled  to  instability,  the  behavior  of  stage  two  appeared  to  contribute  to 


increasing  levels 

of  activity 
Table  8. 

in  stage  1. 

Model  Stage  Initiating 
Model  Results 

Instability 

Simulated 
Reading  No. 

Frequency 

W2/Wi 

AMPF 

Unsteady 

Force 

100%  N//0 

- 

25 

1 

? 

7 

71 

42 

1 

1 

7 

- 

60 

1 

? 

1 

69 

80 

2 

2 

2 

64 

118 

2 

2 

2 

64 

118 

2 

? 

7 

65/66 

118 

2 

? 

2 

41 

220 

2 

2 

7 

43 

350 

? 

? 

7 

80%  N//e 

82 

42 

2 

2 

2 

77 

75 

2 

2 

2 

61 

118 

7 

2 

? 

57 

240 

7 

2 

2 

53 

350 

? 

? 

? 

117 
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At  100%  speed,  rotor  1 operated  with  its  mean  incidence  angle  much 
nearer  steady-state  surge  than  did  rotor  2 and  as  the  fan  approached  insta- 
bility at  low  frequencies,  the  activity  level  of  rotor  1 was  sufficient  to 
initiate  instability  in  the  first  stage.  As  frequency  was  increased,  stage 
two  activity  became  so  pronounced  that  instability  initiation  must  be 
attributed  to  behavior  of  the  second  stage. 

In  the  case  of  80%  speed  operation,  the  mean  rotor  incidence  angle  is 
positioned  nominally  the  same  for  both  rotor  1 and  rotor  2 relative  to  the 
steady-state  surge  location.  For  each  case  that  the  stage  initiating  the 
instability  could  be  identified,  it  was  determined  to  be  stage  2. 

The  stage  initiating  instability  for  the  350  Hz  case  at  both  80%  and 
100%  rotor  speed  could  not  be  identified  by  any  of  the  current  instability 
identification  techniques.  These  instabilities  did  not  appear  to  exhibit 
any  abrupt  or  discontinuous  activity  that  could  be  useful  in  establishing 
the  initiating  stage.  It  is  noted  that  the  results  obtained  from  each  of 
the  three  parameters  given  in  Table  8 are  consistent  with  each  other  and 
tend  to  justify  the  technique  employed  for  identifying  the  stage  where 
instability  initiated. 


5 . 5 EFFECT  OF  UNSTEADY  BLADE  CHARACTERISTICS 

Analysis  of  the  stationary  waveform  data  shown  in  Section  5.2.1  revealed 
that  above  75-80  Hz,  a phase  difference  developed  between  the  model  results 
and  the  test  data.  It  appeared  that  the  PS14  model  results  led  the  test 
data  and  at  118  Hz  the  phase  angle  appeared  to  be  greater  then  100°.  These 
results  seemed  to  be  the  first  indication  that  the  response  of  the  blade 
circulation  could  no  longer  be  assumed  to  be  instantaneous  and  that  the 
finite  time  required  for  the  blade  circulation  to  readjust  must  be  taken 
into  account. 

The  method  of  Section  4.2.3  for  including  the  effects  of  blade  circu- 
lation readjustments  was  included  in  the  model  and  new  simulations  corres- 
ponding to  the  simulations  of  the  Section  5.2.1  100%  corrected  speed,  80  Hz 
and  118  Hz  cases  were  conducted.  The  results  are  shown  in  Figures  63  and 
64  for  just  the  PTIGV  and  PS14  parameters.  The  results  at  80  Hz  show  no 
significant  effect  of  time  required  for  blade  circulation  readjustment,  but 
at  118  Hz  the  blade  circulation  readjustment  introduces  a lag  of  approxi- 
mately 25  degrees.  This  lag  is  not  sufficient  to  cause  a match  between  the 
model  results  and  test  data,  but  the  trend  is  correct.  As  was  stated  in 
Section  4.2.3,  refinement  of  the  parameter  time  constant  coefficient  C3  is 
required . 

Two  interesting  facts  emerge  from  this  analysis:  1)  The  introduction 

of  the  finite  time  required  for  blade  circulation  adjustment  introduces  a 
phase  lag  not  only  in  the  parameters  downstream  of  the  fan,  but  introduces 
noticeable  waveform  changes  in  both  the  upstream  and  downstream  parameters. 
This  finding  supports  the  growing  realization  that  unsteady  flows  in  com- 
pression components  and  associated  ducting  form  highly  coupled  systems  and 
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63.  Circulation  Readjustment, 


100%  N/V§\  80  H/.  • 


implies  that  any  change  to  a compression  component  must  be  carefully  evaluated 
by  examining  the  system  unsteady  response  rather  than  just  an  evaluation  of 
the  local-steady-flow  response.  2)  The  largest  incidence  angle  lag  occurs 
in  stator  2 at  both  80  Hz  and  118  Hz  while  there  is  almost  insignificant 
lag  in  the  other  four  blade  row3.  This  occurs  in  response  to  the  greater 
activity  of  the  rotor  2 incidence  angle  which  leads  to  large  variations  in 
rotor  2 deviation  angle,  and  thereby  directly  affects  the  stator  2 response. 
Examination  of  the  Appendix  C plots  at  80  Hz  and  118  Hz  at  100%  corrected 
speed  will  provide  graphic  illustration  of  this  effect. 


121 


6.0 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  program  has  provided  the  opportunity  to  investigate  some  of  the 
capabilities  of  a dynamic  digital  model  of  a fan  component  and  to  compare 
the  simulation  results  against  test  data.  This  verification  step  has  been 
missing  in  most  of  the  previously  conducted  dynamic  modeling  studies.  Based 
upon  the  results  of  this  program,  it  has  been  concluded  that  the  Dynamic 
Digital  Blade  Row  Compression  Component  Stability  Model  correctly  calculated 
the  manner  in  which  unsteady  flows  interact  with  ducts  and  a fan  component 
and  is  limited  only  by  knowledge  of  the  blade  characteristics  input. 

Further,  it  is  concluded  that  this  type  of  sophisticated  one-dimensional 
model  offers  a unique  method  for  studying  compression  component  stability 
i haracteristics,  especially  where  the  disturbance  is  of  a nearly  planar 
nature.  These  general  conclusions  are  supported  by  the  following  specific 
f indings. 

• The  clean-inlet-flow  fan  map  was  accurately  reproduced  at  80%  and 
100%  corrected  speeds'.  The  prediction  of  aerodynamic  instability 
initiation  occurred  at  the  experimentally  determined  surge  points. 

• Accurate  simulations  in  terms  of  waveform  shape,  amplitude,  and 
phase  angle  were  obtained  with  unsteady  inflow  conditions  at  both 
80%  and  100%  corrected  speeds  over  the  42  Hz  to  80  Hz  frequency 
range. 

• At  frequencies  greater  than  80  Hz,  the' simulations  give  proper 
amplitude  trends,  but  do  not  properly  reproduce  the  harmonic 
frequencies  of  the  fundamental  waveform. 

• As  anticipated,  the  values  of  the  predicted  planar-wave  distortion- 
transfer  pitch-line  coefficients  fall  between  the  hub  and  tip  test 
values  at  both  80%  and  100%  corrected  speeds  over  the  42  Hz  to  350 
Hz  frequency  range.  Further,  the  anomolous  "amplification"  at  80% 
corrected  speed  with  118  Hz  planar  wave  inflow  is  predicted. 

• The  amplitude  transmission  characteristics  results  show  that  even 
though  a fan  component  will  attenuate  planar  wave  amplitudes,  the 
individual  stages  or  the  IGV  blade  row  can  amplify  the  planar  wave 
disturbance  depending  on  corrected  speed  and  planar  wave  frequency. 
Amplification  within  an  IGV  blade  row  is  not  understood,  but  has 
been  noted  in  the  IGV  blade  row  for  simulations  of  other  compression 
components. 

• The  planar  wave  distortion  sensitivities  obtained  by  throttling 

the  model  to  instability  reproduced  the  experimental  data  relatively 
accurately  over  the  42  Hz  to  350  Hz  frequency  range  at  80%  corrected 
speed.  At  100%  corrected  speed  the  sensitivities  are  predicted 
only  for  low  frequencies  and  only  the  trend  is  predicted  at  high 
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frequencies.  This  lack  of  agreement  is  attributed  to  an  input 
data  deficiency,  that  is,  the  relative  total-pressure  loss  coeffi- 
cients and  deviation  angles  in  the  post  steady-state  surge  region 
are  based  upon  extrapolations  which  are  not  directly  verifiable. 

• The  planar  wave  distortion  sensitivity  determination  is  sensitive 
to  compression  component  deterioration  and  thus  indicates  that 
the  clean-inlet-flow  surge  line  should  be  carefully  monitored 
throughout  any  stability  test  program. 

• Accounting  for  the  finite  time  needed  for  the  blade  circulations 

to  readjust  led  to  better  agreement  in  the  phase  of  flow  parameters 
downstream  of  the  fan. 

• It  was  found  that  the  stability  of  the  tangential  blade  force 
component  due  to  the  storage  of  angular  momentum  within  a blade 
row  is  very  sensitive  to  the  fan  aerodynamics.  The  fact  that 

under  some  conditions  this  force  exhibits  high  frequency  oscillations 
as  the  region  of  experimental  instability  is  approached,  may  be 
an  indication  that  this  force  may  be  relatable  to  the  aeromechanical 
behavior  exhibited  by  blade  and  vane  strain  gages. 


• The  aerodynamic  instabilities  encountered  during  model  throttling 
simulations  occurred  where  test  surges  were  encountered  and 
exhibited  characteristics  which  replicated  the  behavior  of  the 
actual  physical  instability. 

As  a point  of  discussion,  it  was  observed  that  as  model  instabilities 
develop,  it  was  often  difficult  to  determine  precisely  which  stage  was 
initiating  the  instability.  At  the  heart  of  this  matter  is  an  age  old 
question,  "Which  comes,  first,  the  chicken  or  the  egg?"  In  the  context  of 
the  initiating  stage  discussions,  the  question  was  whether  the  lightly 
loaded  downstream  stage  (mean  point  operating  well  away  from  steady-state 
incidence  angle  at  surge)  which  was  experiencing  large  oscillations,  was 
the  initiator,  or  was  the  highly  loaded  upstream  stage  (mean  point  operating 
close  to  steady-state  incidence  angle  at  surge)  which  was  experiencing 
small  oscillations  the  initiator.  The  coupling  between  stages  is  quite 
obvious  and  in  this  study  it  was  assumed  that  the  downstream  stage  was  the 
initiator  since  its  dynamics  seems  to  cause  the  upstream,  more  heavily 
loaded  stage  to  enter  and  remain  in  the  post  steady-state  surge  line  region 
for  a significant  amount  of  time.  Clearly,  more  study  and  analysis  of 
instability  development  are  needed  if  the  fan  and  compressor  designers  are 
to  know,  without  ambiguity,  which  stage  requires  refinements  when  stability 
improvements  are  needed. 

As  a result  of  the  simulations  and  the  ensuing  analysis,  it  has  become 
clear  that  the  capabilities  of  the  model  could  be  further  improved  and  that 
a deeper  understanding  of  the  physical  processes  which  take  place  in  a 
compression  component  during  unsteady  flow  could  be  gained.  The  following 
recommendations  address  these  problems. 
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Refinements  to  the  blade-row  time-constant  coefficient/solidity 
relationship  are  required  if  a more  exact  match  between  the  phase 
angle  of  test  and  model  parameters  is  desired. 


• A method  for  extrapolating  the  steady-state  blade  characteristics 
beyond  the  experimental  surge  point  must  be  developed. 

• Examination  of  the  effects  of  free-volume  lengths  and  blade-row- 
volume  lengths  needs  to  be  accomplished  for  frequencies  above  80 
Hz  to  determine  which  is  the  controlling  factor  for  ensuring  that 
the  harmonic  frequencies  are  properly  reproduced. 


APPENDIX  A 


DATA  ENHANCEMENT  OF  CORRECTED  AIRFLOW  AND  PRESSURE  RATIO 


The  results  of  the  data  enhancement  analyses  are  given  In  Figures  65-67 
for  frequencies  of  42,  80,  and  350  Hz  at  100%  corrected  speed  and  In  Figure 
68  for  75  Hz  at  80%  corrected  speed.  It  should  be  noted  that  data  enhancement 
results  are  presented  only  for  those  points  where  the  airflow  measurement  was 
determined  to  be  properly  operating  and  for  those  points  which  correspond  to 
cases  that  were  simulated  (Rdg.  No). 
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APPENDIX  B 


LOSS  COEFFICIENTS,  DEVIATION  ANGLES. 
AND  LIFT  DIRECTION  CORRECTION  ANGLES 
AS  A FUNCTION  OF  INCIDENCE  ANGLE 


The  rotor  loss  coefficients  and  deviation  angles  for  100%  corrected  speed 
are  shown  graphicaily  in  Figures  69-72,  while  the  rotor  loss  coefficients  and 
deviation  angles  for  80%  corrected  speed  are  shown  graphically  in  Figures 
73-76.  The  polynomial  coefficients  used  to  fit  these  curves  are  given  in 
Tables  9 and  10  for  100%  and  80%  corrected  speeds,  respectively.  The  lift 
direction  correction  angles  for  the  IGV,  rotors,  and  stators  are  shown  graph- 
ically in  Figures  77  and  78  for  100%  and  80%  corrected  speeds,  respectively. 
The  coefficients  for  the  linear  curve  fits  of  Figures  77  and  78  are  given  in 
Table  il  for  both  100%  and  80%  corrected  speeds. 
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Stage  Stack  Results 


Figure  70.  Two-Stage  Fan  Rotor  1 Deviation  Angles,  100%  N /J§, 
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Figure  76.  Two-Stage  Fan  Rotor  2 Deviation 


Table  9.  1'olynomial  Representation  of  Two-Stage  Fan  Characteristics,  100%  N/s/S". 
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Table  10.  Polynomial  Representation  of  Two-Stage  Fan 
Characteristics,  80%  N/s/o. 


A)  RELATIVE  TOTAL-PRESSURE  LOSS  COEFFICIENT,  (JU  = P^(TAN(  i) ) 
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TANGENT  OF  INCIDENCE  ANGLE 


APPENDIX  C 


r 


PLOT  DATA  FROM  STATIONARY 


MODEL  OPERATION 


This  appendix  contains  rotor  entrance  flow,  tangent  of  rotor  incidence 
angle,  rotor  relative  total-pressure  loss  coefficient,  tangent  or  rotor 
deviation  angle,  rotor  blade-row  total-pressure  ratio,  rotor  blade  row 
static-pressure  ratio,  and  rotor  blade-row  total-temperature  waveforms  for 
the  case  where  the  simulation  had  established  stationary  operation  prior  to 
the  initiation  of  throttling  to  surge.  The  waveforms  for  25,  42,  60,  80, 
118,  220,  and  350  Hz  at  100%  corrected  speed  are  given  in  Figures  79-127. 
The  waveforms  for  42,  75,  118,  240,  and  350  Hz  at  80%  corrected  speed  are 
given  in  Figures  128-162. 
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Figure  98.  Blade  Row  Static-Pressure  Ratio  Vs. 
lOO%  NATS,  60  Hz. 


tr»/, 


Time 


.04 


.01 


-.02 


-.05 

0 

TIME,  s^c 


Figure  103.  Tangent  of  Deviation  Angle 
Vs.  Time,  100%  N/J&,  80  Hz, 
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Figure  105-  Blade  Row  Static-Pressure  Ratio 
Vs.  Time,  100%  N/v/^,  80  Hz. 
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Figure  108.  Tangent  of  Incidence  Angle  Vs.  Time 
100%  NA/5,  118  Hz. 
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Figure  109-  Loss  Coefficient  Vs. 

100%  N/s/e,  118  Hz. 
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Figure  ill.  Blade  Bow  Total-Pressure  Ratio 
Vs.  Time,  100%  N/V^T,  118  Hz. 
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Figure  112.  Blade  Row  Static-Pressure 
Ratio  Vs.  Time,  100%  N/JG, 
1 18  Hz. 
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Figure  119.  Blade  Row  Static-Pressure  Ratio  Vs.  Time 
100%  N /$,  220  Hz. 
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Figure  120*  Blade  Row  Tota 1 -Tempera ture  Ratio  Vs.  Time 
100%  N/-J&,  220  Hz. 
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Figure  125. 


Ulado  Rov  Total-Pressure  Ratio  Vs.  Time,  100% 
N/-J&,  350  Hz. 
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Figure  13(>.  Tangent  of  Incidence  Angle  Vs.  Time,  80%  N/\/9 
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Figure  147.  Blade  Row  Static-Pressure  Ratio  Vs.  Time 
80%  N/v^,  118  Hz. 
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Figure  150.  Tangent  of  Incidence  Angle  Vs.  Time 
80%  N/v/ff,  240  Hz. 
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Figure  151-  Loss  Coefficient  Vs.  Time,  80%  N/V^ 
240  Hz. 
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Figure  153.  Blade  Row  Total-Pressure  Ratio  Vs. 
30%  N/V0,  240  Hz. 
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Figure  15^-  Blade  How  Static-Pressure  Ratio  Vs.  Time 
80%  N /J§,  2 40  Hz. 


Figure  157-  Tangent  of  Incidence  Angle  Vs.  Time 
80%  NA/0,  350  Hz. 


22'* 


BLADE  ROW  TOTAL-PRESSURE 


l.U 


VWVAAA/W'^/ 


VvVWWV^V 


.01  .02 

TIME,  sec 


.03 


Figure  l6l.  Blade  Row  Static-Pressure  Ratio  Vs. 
Time,  8o%  N/V^,  350  Hz. 
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APPENDIX  D 

AMPLITUDE  TRANSMISSION  CHARACTERISTICS 


The  IGV,  Stage  1,  Stage  2,  and  overall  amplitude  transmission  character- 
istics are  given  for  each  frequency  in  Figures  163—169  for  100%  corrected  speed 
and  in  Figures  170-174  for  80%  corrected  speed. 
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AMPLITUDE  RATIO 


£3  BASELINE 
O NEAR  INSTABILITY 


AMPLITUDE  TRANSFER  LOCATION 


Figure  167  . Amplitude  Transmission  Characteristics,  100%  , 118  Hz. 
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Figure  169.  Amplitude  Transmission  Characteristics,  100%  N/v/57 350  Hz 
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Figure  1 7 1 - Amplitude  Transmission  Characteristics,  80%  N/VS"- , 73  Hz 
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Figure  173*  Amplitude  Transmission  Characteristics,  80%  N/\/0  , £/tO  Hz 
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APPENDIX  E 


INSTABILITY  IDENTICIATION  PARAMETERS 


The  plots  of  the  rotor  ’’unsteady  blade  force,"  amplification  function,  and 
tangent  of  the  incidence  angle  parameters  at  the  time  of  instability  develop- 
ment are  given  in  Figure  175-204  for  all  the  frequencies  simulated  at  both 
100%  and  80%  corrected  speeds. 
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Figure  175.  Unsteady  Blade  Force  Vs.  Time,  100%  N/V^  ^ 
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Figure  170.  Amplification  Function  Vs.  Time,  100%  N/\/§" 
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Figure  179.  Amplification  Function  Vs.  Time,  100%  N/s/0~  60  Hz 
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Figure  l8l  . Unsteady  Blade  Force  Vs.  Time  , 100%  N/\/5"  , 118  Hz 
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Figure  183 
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Tangent  of  Incidence  Angle  Vs.  Time 
100%  NA/e  118  Hz. 
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Figure  184.  Unsteady  Blade  Force  Vs.  Time 
100%  N/v/0  220  Hz. 


Figure  185 • Amplification  Function  Vs.  Time 
tOO%  N 220  Hz. 
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Figure  187.  Unsteady  Blade  Force  Vs.  Time  t 100%  N /-/&  , 350  Hz 
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Figure  188.  Amplification  Function  Vs.  Time  , 100%  N/v/^T , 350  Hz 
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Figure  189.  Tangent  of  Incidence  Angle  Vs.  Time 
100%  N/V0  , 350  Hz. 
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Figure  191.  Amplification  Function  Vs.  Time  , 80%  N/v^T  t 4 


Figure  192.  Tangent  of  Incidence  Angle  Vs.  Time 
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Figure  193.  Unsteady  Blade  Force  Vs.  Time  , 80%  N/\Z5~,  75  Hz 


Figure  199.  Unsteady  Blade  Force  Vs.  Time  , 80%  N /JQ  , 2'tO  Hz 
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Figure  201<  Tangent  of  Incidence  Angle  Vs.  Time 
80%  N/*aT  , 240  Hz. 
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Figure  204 • Tangent  of  Incidence  Angle  Vs.  Time 
8o%  N/n/0  , 350  Hz. 
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GLOSSARY 


Nond imens ional 
Amplitude 


Amplitude  Ratio 


Planar  Wave  Distortion 
Transfer  Coefficient 


Planar  Wave  Distortion 
Sensitivity 


One-half  of  the  difference  between  the  wave  maximum 
and  wave  minimum  divided  by  the  average  value  of  the 
wave. 

Ratio  of  nond imens ional  total-pressure  amplitude  at 
exit  to  the  nondimensional  total-pressure  amplitude 
at  entrance  taken  across  the  IGV-,  Rotor  1,  or 
Rotor  2. 

Ratio  of  the  fan  discharge  nondimensional  total- 
pressure  amplitude  to  the  fan  entrance  nondimensional 
total-pressure  amplitude. 

The  loss  in  surge  pressure  ratio  at  given  frequency 
measured  from  the  clean  surge  line  to  the  mean 
operating  point  at  surge  with  planar  wave  distortion 
divided  by  the  nondimensional  amplitude  of  the  inlet 
total-pressure  fluctuations. 


NOMENCLATURE 


A 

a1B 

A* 

B 

C 

C(k) 

Ci 

CL 

CZ 

cu 

C3 

Di 

FB 

fd 

fdz 

fl 

ft 

IGV 

L 

L 

M 

M 

N 

N//e 


Area 

Area  Normal  to  Flow  Direction 
Choked  Area  at  M = 1 

Y-Axis  Intercept  of  Lift  Direction  Correction  Angle 
Absolute  Velocity 
Theodorsen's  Function 

Coefficients  for  Relative  Total-Pressure  Loss  Coefficient 
Polynomials 

Lift  Coefficient 

Axial  Velocity 

Absolute  Tangential  Velocity 

Time  Constant  Coefficient 

Coefficients  for  Deviation  Angle  Polynomials 
Blade  Force 

Drag  Force  in  Direction  of  Blade  Entrance  Flow 

Axial  Component  of  Drag  Force 

Blade  Lift  Force 

Blade  Tangential  Force 

Inlet  Guide  Vane 

Volume  Length 

Lift  Force 

Mach  Number 

Slope  of  Lift  Direction  Correction  Angle 
Rotor  Speed  in  Revolutions  per  Minute 
Corrected  Speed 


r 


NOMENCLATURE  (Continued) 


Static  Pressure 
Dynamic  Pressure  (PT  - P) 

Mean  Static  Pressure 

Polynomial  of  Degree  n 

Planar  Pressure  Pulse  Generator 

Total  Pressure 

Loss  in  Surge  Pressure  Ratio 

Specific  Gas  Constant 

Region 

Blade  Spacing  in  Definition  of  Solidity 
Surface 

Entropy  Production  Term 
Total  Temperature 
Pitch  Line  Wheel  Speed 
Freestream  Velocity 
Relative  Velocity 
Volume 
Flow  Rate 
Corrected  Flow 

Approximate  Corrected  Airflow  (See  Section  5.2.1) 

Relative  Velocity 

Stagnation  Velocity  of  Sound 

Blade  Chord 

Frequency 
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NOMENCLATURE  (Continued) 


Gravitational  Constant 

Incidence  Angle 

Kinetic  Pressure,  1/2  p (W)^ 

Pitch-line  Radius 
Specific  Entropy 
Time 

Angle  of  Attack 

Blade  Inlet  Air  Angles  (Stator  and  Rotor,  Respectively) 

Blade  Inlet  Metal  Angles  (Stator  and  Rotor,  Respectively) 

Lift  Direction  Vector  Correction  Angle 

Lift  Direction  Vector 

Ratio  of  Specific  Heats  for  Air 

Deviation  Angle 

Small  Quantity 

Efficiency  ('i'n,/^) 

tt/tstd 

Wa ve ] eng  t h 
Dummy  V tr iable 
Density 
Solidity 

Blade  Row  Time  Constant 
Blade  Camber 

C 

*4 

Flow  Coefficient  ( ) 


NOMENCLATURE  (Continued) 


rm 


\b' 


7u' 

U)* 

O) 


Underline  - 


Work  Coefficient  y-1  ^°R 


(TTi+l  - tTi) 


[2wN(r1) J /2 


Pressure  Coefficient  y-1  T^i  [ Ptj 


pt 


(T-1)/y 


i+1 


[2irN(r1)]i/2 

Relative  Total-Pressure  Loss  Coefficient 
Reduced  Frequency  (2irfc/W) 

Circular  Frequency 
Vector  Quantity 
Scalar  Product 
Vector  Product 


-1 


SUBSCRIPTS 

0 

1 

2 

E 

I 

M 

S 

T 

T 

i 

k 

m 

n 


Mean 

Compressor  Entrance  Station;  also  Blade  Row  Entrance  Station 

Compressor  Discharge  Station;  also  Blade  Row  and  Exit  Station 

Effective 

Instantaneous 

Mean 

Standing 

Total 

Traveling 

i-th  Station 

K-th  Volume 

M-th  Stage 

N-th  Blade  Row 
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NOMENCLATURE  (Concluded) 


SUPERSCRIPTS 

Volume  Average 
Relative  Frame  of  Reference 
Derivative  with  Respect  to  Time 
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